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ABSTRACT 


Title of Thesis; Extensions of the Laplace Cascade Method. 
John Hilary Billings, Doctor of Philosophy, 1960. 


Thesis directed by: Research Professor Joaquin B. Diaz. 


The Laplace cascade method is concerned with second order linear 


hyperbolic equations of the form 


(1) Wey t a(x,y) vw, + b(x,y) wy +ec(x,y) u =O. 


The substitutions u,™ uy + au and u_, = a + bu lead to the 


equations 

(2) hy + bu, - hu ™ 0 
and 

(3) -" au, ~kum™ 0, 


where h ™a,+ab-c and k ie ab -c are the two Darboux 
invariants. If either h or k vanishes, (1) has been reduced to a 
system of two first order equations, while if neither vanishes the 
equation may be cascaded in two directions. Solving (2) (or (3) for 
u in terms of u, (or Uy)» and substituting the resulting expression 
into (1) yields an equation for u, (u_y) of the same formas (1) but 
with new coefficients, in general. This process may be iterated, forming, 
a chain of equations, until either the original equation reappears, or one 
of the corresponding invariants vanishes. 

An extension of Volterra's product integral to the non-homogeneous 


system 





n 
a(x) = > ( 244(x)uj(x)) + £,(2), 
(1) ae 


u, (bd) = Yio 3 i= 1,2, oeogit 





is made first. Thon the Laplace method is extended to systems of 


second order hyperbolic equations, of the form 


8) 


Be n 
(5) dy + y iS 3 "3+ val Day J “3+ a re hn i - 
dxdy Jel ja) Sal vig te=1. 


. LH1,2,eceyn e 


Matrix notation is promptly introduced, and equation (5) is rewritten 


as 


(6) u,, + au, + BUY + cu = 0, 


The two related matrix invariants for (6) are 
me A + BA - C 
K = BY + AB = C 3 


and the chain of equations is developed as for the single equation. If 
either H or K is identically zero, the resulting two systems of first 
order equations,may be solved by employing the above-mentioned extension of 
Volterra's product integral. | 

The invariance of H and K are discussed in a manner analogous 
to the invariance of the functions h and k of the single equation. 


This is followed by a consideration of periodic systems, i.e. systems such 





that after j iterations the original equation reappears. This dis- 


cussion results in two theorems, the first of which is 


Theorem I - A system of equations of the form (6) having 
constant matrix coefficients A and B, can be reduced to the form 
8 + af 
“_ = ES 
“Vv 
by a change of variables U = Av’, if an only if AB = BA. The second 
theorem, arising from systems of period two leads to a discussion of the 


form of the solution to the matrix analog of Liouville's equation 


YO _se9_269 
“a” CO 


Discussion of the form of the solution when the chain terminates 
after a finite number of iterations leads to two further theorems which 
are completely analogous to theorems proved by Darboux for the single 
equation. 

The second extension of the Laplace cascade method is made to the 
third order linear hyperbolic equation in three independent variables of 


the form 


(7) yy, * aul, + Bip, b Cl aa, + oR ate +gu=0, 


Here the number of invariant functions to be consider jumps fram two to 
eighteen. The nature of the “invariance” of these functions is different 
from that of the second order invariants, in that some of the functions 

are truo invariants, while others can only be considered as quasi-invariants. 
Four mothods of cascading the oquations are discussed, each of which requires 


severe restrictions on the coefficients a, Dycece,g The class of equations, 





for which each method will result in a termination of the chain after a 
finite number of iterations, is explicitly pointed out. 

The final extension is a generalization of this third order exten-= 
enon. to the nrh order linear hyperbolic equation with n independent 
variables, of the form 


+ 


n 
u : a. 
a Se a ix 


eos) a ae 
(8) ; 
+>. a + ee eee 


yelp yay RH 
i+] 
Introducing the linear operator D, so that (8) becomes 
D(u) = 0; 
the number of identities which may lead to a decomposition of (8) into 
a system of two equations of lesser order is computed. The exact number 
of true invariants is determined to be n(n-1), while upper and lower 
bounds on the number of quasi-invariants are formulated. Theorem V proves 
that an invariant h is atrue invariant, i.e. invariant under the change 
of coordinates uw = \GaieXo2 eax, ) 7 aE 
h= 0 a * SafG) ~ 489 
cos 
for some i,j = 1,2,c+e9M3 i * jy while any other invariant is, in general 
not a true invariant, hence a quasi-invariant. 


A discussion of earlier extensions is included. 
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SECTION I 
INTRODUCTION 


The cascade method was devised by Pierre Simon Laplace, and 
investigated in camplete detail by Gaston Darboux.© Earlier ex- 
tensions were made by U. Dini,3 and J. LeRoux’, and these exten- 
sions are discussed briefly in SECTION VI of this paper. “The method 
itself deals with the linear hyperbolic equation with variable coef- 


ficients which has the form 





(1) fu) = 920 y a(x,y) “du + b(x,y) Ju_+ e(x,y)u = 0, 
| dxaY ax Sy 
where the coefficients a, b, and ¢ are analytic in a certain domain 
D, or at least as many times continuously differentiable as we feel is 
necessary, and u —u(x,y) is-a real function of real variables. 
Leplace's cascade method begins with the introduction of what 
are now called the two Darboux invariants, h and k, which are defined 


by the relations 


Jj 


(2) yb 
k => 92+ ab <-e 
oy : 


h-= 08 4 ab =¢ 
ox 


Then (1) may be written in either of the following forms: 





so (BE sau) + (GE tand-m = 05 


(3) 


ey (oh + Ou, > 2 
53 ; ) a( a bu) = ku = 0. 


We may then consider the obvious substitutions 


al - 22 + gau , 
(4) , 


U7 g U + bu 9 


and further reduce (3) to 


So 
Sal 
om 
— “1 


Now if we should be so fortunate that either function h or k is 


+ bu, » hu = 0 


wo 


(5) 


fo 
& 
i 
2 
" 
© 


identically zero, we will have succeeded in reducing our original 
second order equation to a system of two first order equations, 


either 





(6a) 








‘hs oe 





or 
oo = gu + bu 
Q ) 
6b 
(6b) sea, 
yy on) = 0 ; 


either of which system may be solved by quadratures. 

If, however, as is more likely the case neither h nor k 
is identically zero,all is not lost. To be specific, let us consider 
the system (6a), as the details are “quite similar whichever system 


we choose.” Instead of (6a) » we have 
(6a") 


We may integrate the first equation of (6a') to obtain 


~ \ aay Sady 
(7) u =e ie uy + x(x) 


where X(x) is an arbitrary function of x. Substitution of this ex- 


pression for u into the second equation of (6a') yields 


x 


— \ ady ady 
2 4 bu, - he \. J 3, dy + Xx)|=0, 





or 


Jer fF au 
(8) e + bu - {aay 
he Ox “ \e vu, dy + X(x)- 





Taking the partial derivative of both sides with respect to the vari- 


able y gives 








{aay 
Sy ju 
U 
o.. * t+ =| - d log h O Ll + bul + 
Ox : ay 5 x wv 
3s | 
U 
_ _t tot d i” Ob u 2 aay 


After some straightforward algebraic manipulations, this becomes 


2 
le os. 








ee 1 + iH 
yx oY yx a Fy + cy w= 0; 
(9) 
where a, =a - Olog h b = by, and 
Oy ee 
Cy = Sey oe Oa + Ov ~ b 2 log h 


er: “a7 
We observe that (9) is of the very same form as (1), but with new 


coeffients a by» and Cs Hence we may repeat our original pro- 


lg 
cess, in hopes that either of the new invariants hy or ky will 
perhaps be zero. If not, we may proceed to evolve a new equation for 


the variable ou 93 and in fact we may iterate the entire procedure as 





often as necessary to produce a chain, or cascade of equations, in the 
hope that at some point the iteration will stop because ome of the in- 
variants will become zero. This will enable us to solve a first order 
system by quadratures, and by tracing our way back through the chain, 
we can easily solve the original equation (1). 

As indicated previously, we may do an enslogous procedure with 
system (6b), producing a caseade of equations in the “opposite direc- 
tion". Darboux poimted out that if we followed the h substitution 
with a k substitution, we would not produee a new chain, but would 
indeed revert to the original equation (1). In facet, if we denote (1) 
by E, denote the equaticms obteined from E by use of the functions 
h, hy, ho,» coo by Ej5 Eos EB, coe, and denote those obtained from E by 
use of the functions k, kj, Kp, oo by E -y» 2.2; B 3) cee o, Our chain 
of equations appears as 


EB 


oo oy EB & E, 5); po 2 ° ° 


“2? “al? 


function, and 


If we should take any Bo obtained by use of an boa 


attempt to obtain a new equation using a x function, we would in fact 
produce equation En-1° 
Darboux diseusses these and many others points regarding the 


The same relationships hold for the E . equations. 


Laplace eascade method, including the nature of the invarience of the h 
and k functions, periodicity of the h and k functions, and the form 
of the most general solution obtainable if the chain terminates after a 
finite number of iterations in either direetion. It is the purpose of .f. 
this thesis to extend this eascade method to larger classes of equations 
and to carry out similer investigations regarding these new applications 


of the method. 





In particular we will first show how the cascade method can be 
applied to systems of n linear second order hyperbolic equations in 
two independent variables, with n dependent variables, discussing the 
invariant nature of our substitution functions, the general form of a 
solution when the chain terminates after a finite number of iterations, 


and the solution of some typical systems. We will consider the system 


2 z n 
. u u. 
(10) oe DF a, OCD + "43 os 
0 Xov jai Ox j=1 Oy 
nN 
- 
2 “a5 vs - QO, 1 = Ne fhe ea 3 
get 
a e 
oe "5; - 4:5 &y); bi5 = Dig (%s¥)s Cys = c, (x,y) are continu- 


ously differentiable as often as necessary, and the u, = uj, (xsy) 
are real functions of real variables. We willfirst express this equation 


in matrix form, letting A = (a3 5). Bas (d5 5), C= (cs) ’ 


v=[u]|> (o) = Ons O=[ 00 «46.0 
Us, 0 oe woh OC 
ut QO GO. veicien © 


Then (10) may be written as 


eu + Clee = 
Gd) ae F (0) 











2 

where De. a 3 i U = dy U = ety 
dxay|? 2 eS | oe = 
Q ° a 
>, 2%, a3. 
OxOY OX y 


Introducing the two Darbour nxn matrix invariants 


sz A + BA - C 
¥ Dy 


kK = * + AB = C€ ¥ 
we will be able to proceed with development of chains and solutions 
by the cascade method. 
Our second extension of this method will be to a single third- 
order linesr hyperbolic equation in three independent variables, of 
the form 


(12) Boye + alx,y,z yay + (EF, z yu, + e(x,y,2 yy ss 


+ a(x,y,z2 Ju, + o(x,752 yu, + f(z,7,2 a, + a(x, y,2) u= 


Since we will be dealing in three independent variables, we will see 
that three chains will result, in the "x-direction", the "y-direction", 
and the “z-direction”. Also, sinee the equation is of third order, we 
will see thet sighteen “invariant” functions must be introduced. We 
will note that the invariance of these functions will not be of the same 


nature as that of the eorresponding functions in the second order, two 





must be placed on the coefficients a(xsy5z)5 o o © » &(Xs¥sxX) before 
a chain of equations can be developed. 

Finally we will generalize the results indicated in the preced-= 
ing paragraph to the single linear hyperbolic equation of nt® order in 


n variables, which has the fom 


0 i a a. nel 


0X OXpere 0 x, fie : 


0 X00 9X421 DXi +] 066 OX 


be 


n, ae 
ae 
3 re 0 en a ae 
t3J = if: OX 000% 4 3 Xi ypees 9X, 7 9X5 47006 OX, 
i#3 
Nn 


In this section we will introduce an operator notation which will 

greatly simplify the calculations for ee order equations. This 

will enable us to predict thenumber and form of the identities, and 

the number and form of the invariants corresponding to each identity 

related to the n™ order equation. From this we will be able to in- 

dicate the conditions necessary to produce a cascade of equations, 

which will enable us to tell when a reduction in order will be possible. 
In our concluding section we will discuss briefly the extensions 

made by Darboux, Dini, and LeRoux to systems of second order equations 

in one dependent variable, to a single second order equation in n 


indepent variables, and to a single n*® order equation in two independent 





variables. These extensions are treated more extensively in the works 


of these men noted in the bibliography. 





SECTION II 


AN EXTENSION OF THE PRODUCT INTEGRAL TO 


NON -HOMOGENEOUS SYSTEMS OF THE FIRST ORDER 


In Section III we shall show how the cascade method can be applied 
to a system of second order equations. The end result we hope to obtain 
is a reduction to two first order systems ,one of which will be homogensouw, 
while the other is non-hamogeneous. The solution of the homogeneous 
system ean be found by employing VYolterra's product iseeeaal As yot, 
however, this concept appears not to have been extended to non-homogeneous 
systems. ! It will be necessary for us to do this now, to enable us to 
solve completely the second order systems. 


Consider the system of equations 


n 
oh. * du, s > 
ix Ay, uy + Ff 
a J+ 9 i) 
(1) 
uv, (bd) = Ugg f a 1,25 ovegDy 


where the +A and the f, 


functions of the real variable x, on some non-empty interval of the real 


are given continuous, single-valued, bounded 


line, b& x Sc, andthe u, are n given constants. First let us 
write (1) in the notation of matrices. We denote the row matrices 
U(x)=(u, Uw .. . u,) 
F(x) =(f  f, - + » £) 
u.) 


Uy a (1 6 450 o « 2 2O 


10 





+L 


and the square matrix 


‘= (45), x H ° 


It will be necessary to assume thet the matrix A is non-singular, 
that is that the equations of system (1) are linearly independent. 


In the notation above, system (1) can be written as 


0a) = U(Z)A(z) + Hx), 
(2) ° 
U(b) = U. y 


let P, be any partition of the interval [,¢/, such that 


= Xo << X},< woe SX, = &, and let y be any point in the 


interval /x,,,, x,/. We then define U,, by the relation 


B. ™ U an A( y) (x), a Xyy) ic UL FF ¥ (x, a x.) 
(3) 


mE. AC )(x, -x,.,) + r| + F(} (x -x,_,)- 
Here I is thenxn identity matrix with ones on the diagonal and 
zeros elsewhere. From the form of this recurrence relation (3) we 


readily obtain that 


v= U, | | { aCe | + 


m-1 are m 
a) + D EDL TT [ACR yaa, _pI}} G@- 5) + 
Y=j1 jJzv+l 


+ #8, (xm - xm), 


Cl a 
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Next we consider a sequence of such partitions, {P| such that 
ap m—» 0 , AxX=X,- X_,7* 0. Since all the funetiong con- 


cerned are continuous, we may proceed to the limit: 


Cc 


Cc 
(5) Ue) = lim Uv, = Uy [a4 Fac [a ay, 
b b 


M-e9 
- e C 
In this expression, | A= Vaca )an +1) is the "right" product 
- Dd dD 


e 
integral of Volterra, while F(%) | 4a is the row matrix of 


= 


term by texm Riemann integration 


F(§&) [ a. 


, 
In his discussion of product integration, Schlesinger proved the 


ty Ca) 


of the ‘slements of the row matrix, 


8 
following identity: 


C ¥ 
| A =| As; bd<%<e, 
D > : 
and hence 
(6) | — Po "| a 
' b be 


Putting (6) into (5) we obtain 
c e F ol] ¢ 

(7) Ule) = vo | A + \ F( 8) [a [4 a’ 
b b b/ > 
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Cc 
Since | A is constant with respect to the variable of integration, and 
b 


since the scalar a camumtes with every square matrix, we may factor 


this product integral to the right, and write (7) as 


c ¥ vl C 
(8) U(e) = |UD, + \ F(¥%) | “aH a 

b b b 
This function U which we have derived is a function of the end point, 
c, of the interval {1 db sof o Lf we vary this end point in any interval in 
which the functions a, j and f,; a continuous, single-valued and * 


bounded, then U(c) becomes a function of a real variable which we may 


call x, and hence 


y 


x 
(9) Ux) = |u + \ Fy) | 
db 6b 


«J, x 
a\ de | A 
bd 
We assert that (9) is the solution of (2) and hence, element-wise, the 


solution of the non-homogenous system (1). To see this, we observe first 


= 9 x =. NE] 
that when x = b, | A=] » While { F(§ ) | vray = 0, 
b b b 
x. 
ius Up) = U5: Next we note that each row of the matrix [a is a 
10 . 
solution of the matrix equatim dU =—U(x)A(x) . Therefore 
dx 
x 2 
ad 
ee A sa : 
x | | A) *@). 
b b 


Thus differentiation of (9) yields 





Lh 


x ¥ x 
aa = |v 4 ) F(¥) | Nea] A + 


b b 
x 3 1 me 
+ \U. ( (4) | A) ae je a Ap = 

b 0 ax ‘ 
b 

x | = x ¥ -1 x 

= ¥(2)(| A | A + |U, { wR) Al\a& A\A= 
b b ; b b 


= F(x) + U(x)A(x), 


which proves that (9) is the solution of (2). 
Schlesinger showed 11 thet if A is a matrix of constant 
elements, then 


x 
i. {Ady 
[a= e> ; 
b 


Therefore in the special case when (1) is a system with constant 
coefficients, (9) takes the more familiar form 
' 
~\ Ad Ad 
(aay (Ady 
(10) u(x) = | US + \ F(¥)e a |e 
b 


3 


and this’ sointion can be verified even more readily by direct dif- 
ferentiation. We note the analogy between (12) and the solution of 


the single non-hamogeneous first order equation 





i 


du 


ax OT Ox) u(x) + P(x) 


u(b) = a 


This equation has the familiar solution 


x 
J a(r) ay x _(*acgd 
= : t 
>) Ra (x, + \ e Ny £(¥) a¥ ), 
b 


a form which we could obtain fram (12) by taking the transpose, 





SECTION IIIT 


SYSTEMS OF SECOND ORDER LINEAR 


HYPERBOLIC EQUATIONS 


A. Consider the following system of second order linear hyperbolic 





equations: 
(1) yu a ) - 
4 + u u 
Yx0y — 44.2 ¢ 25,25 * 





Cute 
[ 
t- 
fete 
Cur 
QV 
Cur 
it 
a 
f-to 
Cs 
Q 
| 


+ 


nu 
+ a = 
¥ fig se (OP ee eee, 


where the coefficients Bs 4) Dy 5) and Ci; are all real functions of 
the real variables x and y, continuously differentiable in both 
variables as often as necessary, and u,= u,(x,y) are real functions 


of the real variables x and y. To put this system in matrix notation 


3 


let A=(a;;), B=(by;), C =(cey3), Us["2), (0) = |°|, ana 
“3 : 
5 
Up 

we obtain 

2 U = 

(2) = “- ADL. r BU, + cu (0) | 


16 





ae 
ye 
xQ 


a 
ant 


matrices A, B, and C are all non-singular. 


be 


Where “ZL = See = 
axa Y 


, etc. We will assume that the 





Le 
eq 


We introduce the two Darboux matrix inveriants 


(3) 
K = B + ~ AB oC. 


Then we note that (2) can be written in either of the forms 


o fdu + Av) + 2 ( 


ox Oy, 


Q/ 
CG 


a av ) - HW = (0), 


= 


(4) 


A es i: Bu) + b (2g opera a, 
Sy Noe Ox 
If we consider the substitutions 
U = Jv + 
1 y a 
Gog ue c BU ; then (4) becomes 
uv, + BU, - HU = (0); 
x 
=) 
U, + ew (0). 
y 
_foo <¢ 
If now either H = 0 or K = O, where 0 =], 64... our 
Ha ; 
00 Olnxn 
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original system (2) will be reduced to two systems of first order 
equations, which can be solved by the product integrals of Section II. 
We will make the assumption that none of our coefficient matrices are 
Singular, except the matrix 0O. 

Suppose that H =O. Then we have the following systems to 
solve ; 


U + BU = (0) , 


(6) 
U = vu + AU 


The first equation of (6) has the solution 


¥ =] 


(7) a [BC ¥ ow)aF Yr), 


where Y° (y) is a colum matrix of arbitrary functions of y, and 
—1 
( [a 8 ; y) d *] is the inverse of the Volterra product integral 


of B. Knowing U we may then proceed to integrate the second equation 


ar? 
of (6), to obtain the solution 
= a 
(8) U(x,y) = [ale 9 1) 3° (x) + 
Yr % x = 

+ [ac-ray( [ace ¥a9) Yes 
where “ (x) is a column matrix of arbitrary functions of x. 

if 4H = O, but K = Q, we can solve the system 
Ugg “SACI es (O) , 


(9) J 
Uy, = Ue 





1 


by the same methods utilized to solve system (6), to obtain the 


solution 


hae (T = zit) ‘| #0 + eek wh 


+ { a ry) a4 (| ‘sina, x"(5) 4] 


Consider the system of three equations in three unknowns 


Illustrative example No. 13 


Jy 7 du, , ye" Me4F IY; 4 9,47 +e Uy ape uz = 0 | 
dxOy Ox Sze * OF 


(a) y@ Up- © ‘ag. tat 35% Ye - Yup seu, =O 
ws 2 oT Dea Dy “sane ol a po =o. 


This system can be written in the form (2) with matrix coefficients 


SS y ye” y » BS =Bil o o| , Saas 
-e y ey o. lio 
2 
J e” ye" So Oar 


ae 


If we compute the invariant H, we find H = 0. Thus we may reduce 


to the two first order systems 


() he oe eS 


The first system of (b) has the immediate integral 
“xX 
a ae £@) 
where Y¥(y) is a colwm matrix of three arbitrary functions of y. 
Using this expression for U, in the second system of (b), we obtain 


the solution 
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U(x,y) = ({se 1a a | X(x) (heper) ec Hay 


where X(x) is a colum matrix of arbitrary functions of x. 


B. Consider now successively the substitutions 


Y= 8, 
= 0G) 
y = Wy") j 


; : y 
=m x' 
7 ) 
where A is an n x n matrix of twice continuously differentiable 
, 
| / U 
functions of x and y, such that |A|# 9  , and v= | 2? 


e 3 


ay 

n : 

A ug (x,y). Under the change of variables U = Nv (2) becomes. 
A°xy + eax N y)0,, + (BA +* Ax)U, is 

+ (CA + AA, + BA * \ zy) 7 = (0), 


: al 
Multiplying from the left by /\ we obtain 


Br + (AOA 4- RAs. + (A BA +N A oY de 
+ (AC cA Ata Kier Ne: iS + = xy )U' = (0), 


which is of the same formas (2): 


t ‘Trt <= 
(12) Us + A'Ut, + BUY, + Ger (0), 
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where A' = A= aAlh\+ Ni hy j 


(12) Be ONG BATA 
c= ReAt Nan. + A- TBAy tN Ay 


Let us compute the velues of the corresponding Derboux invariants 


H' and K’ for (11). 


‘= ' 1 ¢ vas 
H At 6B b= CaS 


=(N TAK + NN) + (AMBALA A ANtN * Nege 
~(At%eA+ AMa Nt NBA + NAY), 


NT EN 
2 (AA = AIM 4 BAAS 2, ogg 
(33) ox Ox ox oa — 
Wea - AW IAN, 
ax ax 


Using this identity13 we may differentiate the first term of me 


to obtain 


t 


B= -A*A_ A ‘A+ NM Ne NAN, ai WN Ag t 
AAgy tA” BAA+ A By  TAATAA +A oe 
"AiAgy7 Nich - NAA, - ABA, > Ae 

= N'(K, +BA- c\A 


NHN, 


li 


In similar manner we may compute the value of 


t t 


“1 
Ko BY +t AB =a A Kah 
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We see then that this change of variables U = f\ U' merely produces 
@ similarity transformation on the matrices H and K, and indeed, if 
/\ is chosen so that A H=H/\ , then H' = 4, and shoupd 
AK = KA , then K' = XK, 
If we make the change of coordinates x =Q(x), y= \p (y'), 
then the resulting equation 
Osrye + Y (Pave +P (xpBt ee O(a) YW (hou = (0). 


also has the form of (2). Here 


s= Ya 
B' = P(x) 5 
ch = Pl (x) YW’ (se, 
and the correspanding matrix invariants are 
a= O° (x) YX &, 
c= @ (xh W’ (yh x, 
Finally the change of eoordinates x= y', y= x’ merely has 
the effect of interchanging the invariants, so that H'= K,K = H. 
Thus We see that the term “invariant” is correctly chosen. 
C. In the likely event that neither H nor K is 0, let us consider oe 
the system : 
Uy x Uy + AU , 


(14) 
“> BU, HU 


Solving the second system of (14) for U, we obtain 


i wel “1 
(15) : uU = EY Uy, + EB, 


Utilizing tdsentity (13) we further obtain 
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(16) Uy= -—" BE Uy, + ye Viyy - q™ HA BU, zs H” By U1, +H BU, 


Then mutliplying the first system of (14) from the left by H, and 


substituting (15) and (16), collecting terms we have 


+ #70, + BLU + €,8, = 19), 


1 1 ly 
where 
=] “1 
A, = (HR — H,E  ), 
(17) By = 8B, 
I . 
GC, = (HAN - HB )B +B, -H = AjB, +B eg, 


This equation is again of the form of (2) and we may itérate our 
process in an attempt to reduce this system to two systems of first order. 


Computing the matrix invarients for system (17) we have 


Hy ™ Aly + BL Ay CC) 
wo a «lL oJ] 0 ‘ 
m= (HAH ~ HE), + B(HAH - Bu) = (HAR - BHD - 
- By + F 
wn . el «] «al, 
= OH -~K «A, + AB - BA + B(HAH™ oi 8 ) - (HAH “~ HE B+ 


=1 -] 
+ (HAH )g oa (He )y} 


i  E, 
In a similar manner we may solve the system 
we i BU , 


Ujy + AU, = HD, 





eh 


to derivs the system 


(19) Oe © AU + a8, +c_u, = (0) 


where now 


Ay = A, 


Ss ol -1 


CAs (eK - KK" )A -K+ A= 


= BlyAy 5 ‘i. Ke. 


S 


For (19) our matrix invariants will be 


Hi = Ky 


= =“) “1 1 “1 
Kk, = Kk -H ~ By, + BA “AB + A(KBK - KK) - (KEK -KK JA + 


+ (KBR) -(QK™)y 


If we designate the original system of equations (2) by &, 
we have been able to transform E into a system of similar form, which 
we shall designate Ei, by the change of variables UF, = U_. + AU, and 


Jy 
into another auch system, to be dssignated by Blip by the change of 


variables Vly = U, + BU. Provided that none of the resulting H or K 
invariants is non-singular, we may continue this procedure in both 
directions, establishing a chain of systems 

°o ec @ Biot 1p E, ByE,, o ee °0 S 


The question naturally arises, if after making the substitution 





2) 
U, = Uy + AU and discovering that H # O , we transform system 
(2) into system (17), what results if we then make the substitution 


4 = Uy + BU, ? . Presumably we shall again obtein a system of 
the form of (2), but what exactly are the coefficients of this new 
system in relation to the original system (2) 2? Utilizing the expres- 


sions for the coefficients of (19) and recalling that K,=4, we 


obtain 


. — A, = HAE ~- HH 
Leal 1 y ) 


ue ol -1 


= 5B qH = 


rviat AL - 


= H(A, + BA - BEE, ~ a aE” + 
1 1 a mT 
+ Heat {HE HE -BUE 4g, 


Now consider the substitution vu, .= Au’, ana for /\ let us take J, 


lel 


leew. U = HU, From (12) we see the resulting system, again of the 


lel 
form (2), has the coefficients 


q 
= -1 = = 
A= ata b+ EUR = A, 
8 ow) ooh = 
acl HY 3 


c= g7 c, 43 t+ = A, Hy + xe B, Hy t a” Hy = ©, 


Thus we have transformed into a system which is exactly equivalent to system 


( 2); since 
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H= y= pon og 
*, 4 ee) 
7 
HHH Hoo 
| ) 
A te, 
eee - 
| ence a 


This result may be immediately generalized as follows: If, in system 
E, we make the substitution Uy a1 ay +A,Uy, we will obtain system 
Byy4,, but if we make ths substitution U,; i= Uy * B,U;, we will obtain 


the system E for all 1 = secoy +2, -1; 0, 1, 2yee. . The foregoing 


{=}? 
tacitly assumes, of course, that the corresponding H, and K, are neither 
equal to O. 

Remark fl: We may define an exponential matrix eX, where X is 


any nxn matrix, by the formula 


x 2 of 
Bi Salt tse cane + 22. 
2 36 izil i§ 


Now unless the matrix X cammutes with its x - derivative, we cannot 


assert that ade o eX ax , since 
Ax ax 
a ‘ 
QC US THR (E+ Fe L (Mey yRr+ Pax) ,,, 
dx ax ax 3! a& dx ax 


But if X=Ay, say, where A is an nxn matrix of constants, 


and y is a sealar matrix, then dX =A dy , hence 
dx dx 


XAK = AyAdy = Ady « Ay=@Xy, and thus 4 (e*) = e*%& 
ax dx ax dac ax 


Then consider the special case of system (6), when the matrices A and 


B are both matrices of constants. An integrating factor for the first 


(Bax_ Bx 


system of (6) would be = @. Then we may write 
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H= H = HY SE, 
-| — 
eS ey 
) 
ham, 
5 Ri eer - h 
, ence i 


This result may be immediately generalized as follows; If, in system 
E; we make the substitution U; a1 Vay TAU, we will obtain system 
Ex+), but Jf we make the substitution Uy _,= Uy, +B,0;, we will obtain 


the system E For all 1 = sess ~2, -1, 0, 1, 2reee - The Faregoing 


j=)? 
tacitly assumes, of course, that the corresponding H, and K, are neither 
equal to OQ. 

Remark fi: We may define an exponential matrix e%, where X is 


any nxn matrix, by the formula 


Xx 2 i) 
i Tee et le te pas el + De. 
a 38 i2zl i! 


Now unless the matrix X camutes with its x - derivative, we cannot 


assert that d © Me ees dx , since 
an = 


xX 
= = TeE (RX + Hak) + 2 (MY 
dx ax 3! dx 


X+ x-dx) 


Ri 
Rls 
RI 


But if X=Ay, say, where A is an nxn matrix of constants, 


and y is a sealar matrix, then dX =A dy , hence 
dx dx 


LAX = AyAdy = Ady « Ay= 4X y, and thus 4 (e%) = e*ax 
ax dx 4x dic dx 


Then consider the special case of system (6), when the matrices A and 


B are both matrices of constants. An integrating factor for the first 


(Bax_ Bx 


system of (6) would be = @. Then we may write 





eT 


Bx Bx 
U =e 
€ Uy + © B= (0), 


a ¢p>* e Bx 
a ( ~ U,) (0) ) $9 e Uy = ¥(y), 


where Y(y) is a colum matrix of arbitrary functions of y. 


Substituting vs al 4 ad y (y) into the second equation of (6) gives 


~BX 
&e Yy) = 0, + Av. 
(aay _ Ay 
An integrating factor for this system is then € => € 0 


Integrating this system gives the solution 


y= oh | fe” =- Bx tye alt 


eg Mint Bx le "T v(y)ay + x(z)| 


where X(x) is & column matrix of arbitrary functions of x, 
Now consider system (14); where no longer is H= O , and 
assume ones more that A and B are constant matrices. We my inte- 


grate the first system of (14) in the above manner to obtain 


J = en | fe ne + x(2)} 


Substituting this expression for U in the second system of (14 ) 


leads to the equation 


u, + BU, = 5 ew | fe” Uydy + x2)], 


Multiplying fran the left by e 7 x 7 


» this becomes 


“el Ay 
e iat BU, ) . U,dy + X(x). 
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If we now differentiate both sides with respect to y we obtain 


Ay _-l Ay -l 1 
(@™ ap -@™ Bae )(u,_ + BU) + 


: = @ AY 
ARE at BU, + BU) =e™ Uy . 
Multiplying fron the left by ey » and collecting terms we 


finally obtein 


By + AL Vay + BL VE + Cy Vy = (0) 


lay 


Where the coefficients Aj, By , and Cjare exactly as obtained 
for (17). . 

The curious thing about this second method for starting a 
chain of equations is that even if A and B do not conmute with 
their derivatives, we may still use the integrating factor evAdy and 


ao HAV Ray 
dy 


operating formally as above, pretending that the equation =e ° 
Were actually true, we will arrive once more at system (17) with the 
correct coefficients A, , By, and C, . Thus we obtain a valid 
equation from en invalid operation, but an operation which appears 


valid on the surface, and is campletely analogous to the operations per- 


formed on a single equation of this type. 


D. If the chain has been continued to the i+ 18% system of equations, 
we can readily establish the following recurrence relatioabips for the 


matrix invariants: 


= 


ay -1 
(20) Hag 1 Shy > Ky ~ Ay + AgBy ~ ByAg + ByCHpAyHs - Ey iy yaa 


~(HjA;H, ~ - Hy Hy )B, +(HyA,H, “)y - (HiyH Dy 
y ) 
1+1 5 ) 


pL 5 ©o°0 HO » «1, O, uy Cyeee E 





eS 


These may be solved for H, and K, to give 


H,=K 

(21) a L+4¢1 j ; 7 
a —— Ku nae 
Ky= 2K 4. ~ By yy 7 Aah AsBy-ByAgt BilEgyy Akay “Ate ae) 
ae = oly) =] 
- (Ky) Kaye Keer Baer Bat (Kay 4aBay7 (Kaya at1 dx, 
, ae, 
i oe Cs nC, ol, Q, ls Oy obo 

From these we immediately cbhtain the relaticns 


oo} “ay 
aye er’ Maa” “Pa ~ Aaz* As®2 = BA, +B (HAH, - Hy H, ) - 


a =. a -1 
-(HyAyH, Hayy Wey t(HyAghs ), - Gs By Dy 5 
(23) m es <i a 
Byer" By tH + 75, /AjB) - ByAy ~ Apt BAHAY - Hy HS) - 


ay 


x 


oa: allows el 3 
HAH) - Hy Hy UB; HHA sy - (Ah By 


1 = 5 a6 = 2, aL, 0, ky Cy te 


Fimally (5) may now be written more generally as 


U, + BU, =H U 


i 4 dice Lat 
(2h) x scicaaes> 


U + AU K ee U 
~1y <a © a(ja]) ~“-(iel). 
Considering the first system of ( 24+), suppose B is a matrix 
of constant elements. Multiplying both sides by the exponential matrix 


B 
e t (see remark 1.) and solving for Us _j» we obtain 
= ol “Bx Bx 
O44 = Hea C ma (S Us) é 


This immediately leads to an expression for U in terms of U, and its 


derivatives with respect to x, namely 


el  =Bx Bx «lL =B i Bx 
(25) U=H €© os (e Hy e "'e (e"n, oN obo 2E U; Me»), 
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(Bax (Bax 

In the case where the matrix B is such that e (e )= Be 
does not hold,the solution for U in terms of U, becomes much ~ 
more involved. Solving the first equation of (24) for U,_,, we 

a=. “ -1 . 
have U,_) = By iy, + BU, ). Then Uy.o = H,_o (Wana + BU, _4)= 

@1 ) -1 =1 

pic 2 (2 (Hy1(Uy,+ BU) + BH, * (Uy + BU,)). Iterating thic 
process we see that eventually we will obtain U,.,; =U in terms of 


U,; and the i invariants H, H,,...H Using the relation 


i-1° 
Ya Vie, © Sra Beal 

we may, in a similar manner, obtain an expression for Uy, in terms of U. 
E. We shall now give consideration to what may occur in the form of the 
invariants as we iterate along the chain of systems, Perhaps the most 
natural point to investigate is that of “periodic” systems, i.e. systems 
such that after j iterations, we obtain E 3 =E. We shall say that such 
systems have period j. Consider first a system of period 1. Then. 


E,= 5, H, =H, and K, = K. Hence from (20) we have H = K, and 
H=on -H-A +AB-BA + B(HAH -~H\B™)-(HAH” = HyH” )Bt 
on x yZ yi 
. 
=| =. 
+ (HAH ). ~ (HE ‘e. 


In the special case where AB = BA and B(HAH™ - HyH™)=(HAH-l-uu"1)p, 


this reduces to 


(HaH“)), = (Hy), - A, = 0, 


a 


Which has the immediate integral 
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(26) HAH™* - HYH™> - A= -¥(y), 


where Y(y) is an nxn matrix of arbitrary functions of ¥. 


If in addition HA == AH, (26) is further reduced to’ 
@-1 
HH = x(y) 
= 
Hy Ase ; 


This may be further integrated by the product integral to give 


(27) a= ([ Yap) X00 : 


where X(x) is an nxn matrix of arbitrary functions of x. If we 
select for Y(y) @ matrix of constant elements, D this solution 


can be written 
(28) H= 67% x(x). 


Suppose we have a system of period one and let us make a change 
of variables U= f(A, U'. This will not change the periodicity or the 


equality of the invariants, for, as we have seen, if H = K, then | 
t ol ~ J, 
Heme AD HARFN xcAe x 
Let us select /\ however, so that our coefficient A'™0 . This 


simply requires that (/\ go me A. , and hence 


(29) \ = [ ( Alay x(x) , 


where X™ (x) is a colum matrix of arbitrary functions of x. 





se 


t ’ : f ' “-~ . oe ' 1s i? a et :- ' 
Thn H=~A,+ BA - C= ~-C = K, but EBL +AB-~C EBL -C, 
hence ‘By = O . If we make further assumptions regarding the 
character of A and B, our system can be reduced further. Let A and 


B be matrices of constant elements. Then (29) can be expressed as 


~Ay 
(30) N=e- x* (x), 


Since B= AB A + RAs » We may also attempt. to fina /\ such y 


that B = 0. This gives 


“PX 
- IN =e iy) 


Comparing (30) and (31) and noting that e'e@®% = e@P +2 i¢ ana only if 
PQ = QP, Ut we may take for our arbitrary matrices x* = Pe F y = ety 


to get the matrix 


(32) nye * F8 


if and only if AB = BA. 


Using (32) to change coordinates under U = Av » our resulting coef- 


ficients are \. 


and our system (11) has the reduced form 


t tf 


(33) U = HU 
xy 


which is the matrix analog to the telegraph equations. Since, if A and 


»] 


B are constants such that BA = AB, it follows that H = K, we have 


proved the following: 
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Theorem I. A system of equations of the form (2) having constant 


matrix coefficients A and B can be reduced to the form U..= H'U' 





nge of variables U_ 


Illustrative example #2: 


as MY ou, 34 29°34. T1423 ,28 ou, +10u,+45uptligu,= 0 | 
aeoy Ox DE Oe OOF CFCS 


s - 4°04 - 2Mo 3 > Spay 4 Lguy 4129uz; +20u, - 1Tu,g +h9u, = 5 
Sees as oF. FT OF ) 


y%, ,o">4 WY 4 29,4 20up 4 Lidvs + 6u, + l9ugt 79, = 0, 


dxDyY Ox se oy oy ay 


This system may be written in the form (2), with the constant 


matrix coefficients 


A=[1 3 5 B=l7 2. 28) C=—Ti0 4 149 
2 = 1 =a 17 12 20 «LT 49 
oo.) «OT 2 2 sag 6 19 oe 


We observe first that det A = -23, det B =529, and det C =-8,792, so 
that all matrices concerned are non-singular. (Notice that B= A= and 


C=A- - I.) Computing the values of H and K we obtain 
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ao Fe BA - CH I000 Jl 8645 «(149 1O 5 149) 
0001+ 120 -16 kOi~m {20 -17 49) 
000 6 19 80 6 19 79 
Oe = it 
a 10 16 
O°O'™# 
K=By+AB-C = 000 Jl 45 1k9) | jlo 645 Oho] 
000{ +l20 -16 49 20 <17 49) 
000 6 19 80 6 19 79 
ap (ONO = = ae 
O10 
ool 


Thus H= K =I, and we further observe that AB wBA. The conditions 


of theorem I being satisfied, (c) reduces to 


(4) = fv, 


Under the transformation U = Av' » however, we know that 

ai -1 “1 
p= A rh= \ IA=A Acts, so that (d) is in reality a 
system of three “uncoupled” tele, -aph equations, 


(e) iy =U. 


Each equation of this system has the sams solution u', which can 


I 


Having determined U » We must then compute the solution U, using 


Ss 


jw} 


i 
be found in the literature. ? Thus (e) has the solution U = 


the relation U=/\U'. Equation (32) tells us that 


/\ 


“(Ay +Bx) oO Lh 5 2 TF 
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hence our solution is 
-(Ax + By) 
UU = 6 U;, 


F. We consider next the systems of period two, so that Eo = E, 


and hence Hp = H and Kn = K. Equations (20) and (22) tell 


us that 
(34) 2H = 2K = Ayy +A;By - ByAy + By(KA,K™ - KK) - 
Eh é s é 
“(KAK™ ~ KK )B) + (KAXK™), ~ (K,K™), 
2 
and 
ek = 2H - A, + AB- BA + B(HAE™ - HEM) - 
(35) 


- (EAH™ = HyH™)B + (wan™+),- (ayH™),. 


Summing these two equations we obtain 


feu, ape foe, =a) + (aaa = aalyfomage 
al +} aK KK) = (KAK™ - Kyk™)By| + { B(HAH™ - HYH™) : 
(ua = matey = OH, aH, =O. 


Looking at the terms within each pair of curly brackets, we note that 


under certain obvious conditions of commutativity this system reduces to 


(37) (KK). + (HET), = O, 


Equation (37) has the inmediate integral 


(38) Kk + BE = vy) 
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where Y(y) is a square matrix of functions of y. 
If we impose another condition of cammtativity, whereby 
KH = HK, KH = EK. then an integrating factor will be miltiplica- 


tion fran the right by KH. For this leads to 


fe + HyK = Y(HK) , 


2 (eK) = Y(BK) ,” 
a7 


Hence 


(39) mK T (xy 1 | xx), 


where X(x) is a square matrix of functions of x. 
These funetions of x and y are in fact determined by the coeffi- 
cients A, B, and C. Thus we have proved Theorem I: 


if @ sy ations of the form (2) 






following cummtativity satisfied: 


(a) KA = A,K 

(b) HA = AH 

(c) AB = BA 

(a) A,B, = BA 

(e) Bk = KKB) 
() wp" <a> 


(g) Ki = HK 
(a) KE = m, 
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mK = y'(y) x (z), 
where Y'(y) and X'(x) are square matrices of functions of y and 
x respectively. 
Since H is non-singular, there exists a matrix ©) defined by 


16 Now if the conditions 


the relation H= & » so that © = log H. 
of theorem II are satisfied; and in addition, (39) takes the simplified 
form HK = I, then (34) can be written 


2g 2& - (KK), , 


or 
| = 2 
2% - 2F = as 

hence 
(10) re” ae = an), 
Furthermore, if © has the property that © WOW - BOG, 

© Oy Oy 
then H, = e@ J@=H OO , «0 that 
ay y 
(41) 09 — y-lg 
Oy y° 


Substituting (41) into (40) yields the equation 
& ~ 
(42) 2*@ = 2¢ 2€ o. 
OxOY 


and this non-linear system of equations will then have the solution 


© = log H. Equation (42) is the analog of Liouville's 
“Li 


equation. 
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G. We come now to a consideration of the general form of the solution 

when same H, =O H, + 0. 0O€j<i . That is, we iterate the process 
forming the chain of equations E, Ej, Ep, ... B;, until H,;=0O . It 
Will then be impossible to form equation Ey ty henes E, is the first 
system for which the H invariant vanishes. The first equation of (4) 


then has the form 
pu Ay v4 By ( Wy + Ay0,) = (0), 
Oy ay 
This has the immediate first integral 


=, 
(sat a)-([ ete)" an, 


J 
where Y(y) is a colum matrix of functions of y anly. Upon inte- 


grating a second time, we obtain the solution 


—~— 


Hyco Ay (x, ” ) ©) fe +: ae 
| +f { [aster fT 3,838) 1) ay 


where X(x) is a colum matrix of functions of x oly. 


(43) 


This solution is of the form 


w= Q |x(x) + ("6 Y(v) asf 


where QA. and (RQ are n x n matrices of certain well defined functions of 
x and jy. 
Using the iteration scheme of paragraph D, we see that the solu- 


tion for U has the: form 


(4s ) u= ¢ (x + fe vdy) + a, ( r+ \28ray) + 
te Q(x" +\ 26 vdy) + at (xs {xe ry), 
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Hered QA, , eo A, designate certain well defined nxn matrices of 
functions of x and y, and are camposed of certain products of 
H," ,Qj, B, and their derivatives. We note that the colum matrix 
Y always appears under the integral sign, in the most general case. 
If the boundary conditions are such that Y = (©) we see that the 


solution is of the form 


(45) uU= Ax + Aix = Q.x'+ cont a,x), Q, #9, 


and this is the most general solution in which appear arbitrary functions 
of x, and no integral signa. 

Conversely, if the original system of equations hes a particular 
solution of the form of (45), then the repeated application of the ces- 
cade method will certainly lead, after a number of iterations not greater 
than i, toa system for which H, =O. We will now prove this, i.e. 
if (45) is a particular solution, the number of iterations necessary 
before any H-invariant vanishes,is not greater than i. If we substitute 


(45) into the system (2), we will obtain an expression of the form 


‘ xyiit1) L 
Wit WMxre...t MM, ox = (0). 
From the arbitrary character of X, we assert that 
(46) #, =O j =OMenr., 1+i 
If we campute the form of }¥4, and ¥; 41 We obtain 


Hs Gi, i , IAs 
i = t AAs a A Sot = Sa 


hia a AQ, 





LO 


Equation (46) implies that 


(47) Oi = 
Sota Ge, = 6, 
thus M,;= 005.) 4 AAs iy 2 (2s + Aa,)+ (Be r Aas) . 
> F dz Vy oy 


~(4, + BA-C)Q, = G45 4 44,5 - HAL= 0, 
oy 
Hence we have that 


(48) 2Q, 4 + AQ, = FQ, 
Y 
Recall now the first substitution 


(49) a 
ov : 
If we substitute the expression for U given by (45) into (49), 
and whilize, (47), we see immediately that the resulting expression for 
Uy «will have no derivatives of order greater than 1-1. In addition, 


4. 
(48) informs us that the coefficient of x! 1) in the expression for Uy 


will be HA 4° Whenever H vanishes, therefore, the coefficient of 

x (1-1) in this expression vanishes also. That is to say, if H vanishes 
the expression for U;, would have no derivatives of X of order greater 
than i-2. As a consequence, repeated application may lead to a system 
Es for j < i, for which the invarient E, =O ; otherwise we will 


obtain an equation E,, which has a particular solution of the form 
Uy = Ax e 


Yo summarize, we have proved the theorem III: 





Ky 





i # O then the method of Laplace will lead, after a number of 








H. Consider now exprsssions of the form 


(1) 


} 
U = Ax+QX ot. eoe + A,% 3 


which contain a matrix of arbitrary functions of x, and derivatives 
of this matrix up to a specified drder. It is obviously quite often 
possible to express U in terms of derivatives of an order greater than 


specified. That is, if X is expressible as a sum of matrices: 


x= 6% + CX + ...4+KX Y), 


where 6, Oye + ® are certain matrix functions of x, and X3 is a new 
matrix of arbitrary functions of xf then U may contain derivatives 
of X, up to order it KM » Conversely, it may be possible to reduce the 


order of the highest derivative appearing in U. For instance, if 


us Q(x x) + B(x'+ x), 
then the substitution X, = X + X will reduce the order of the 
highest derivative by 1. 
We will say that the matrix U has rank i+1 with respect to 
x if the highest derivative of X appearing in the expression for U is 
of order i, and it is impossible to reduce this order by substitutions 


of the type described above. 
We assert, that if the system (E£,) is the first for which the 





— 
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invariant H, vanishes, then there exists a particular solution of 
the original system, composed of a colum matrix X of arbitrary func- 
tions of x, and derivatives of this matrix up to and including order i. 
This solution is irreducible in the order of derivatives, and hence has 
the rank i+¢1. 

For if there were 4 substitution which would reduce the order 
of the highest derivatives to 1 -/A), say, then there would be a 
system (E,) j=1i-yW < 1, for which this new expression would be a 
perticular integral. But this would imply that H, = 0, which is con- 
trary to our hypothesis. As a result, it is evident that if (E, ) is 
the first system for which H, = 0, then no other system (E ie 7” with 
positive or negative indices will admit a solution of rank i+1 with 
respect to xX. 

The results of the preceding discussion apply without modifica- 
tion to the second substitution, i.e. the substitution which results in 
the chain (E_}) (E_,) eee If this chain eventually terminates at some 


system (E_ 5) for which K_, = 0, then there exists a particular integral 


J 
of rank j +1 with respect to Vy 


U= BY +Q,y + ...t8, ¥ (3) 
and conversely. 


I. We are now in a position to construct all the systems of a given 
dimension{where by the dimension. is meant the number of dependent vari- 
ables) of the form ( 2), which will lead to a general integral by this 


extension of Laplace's cascade method. Suppose, for example, we wish to 
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eonstruct a system of dimension n, which will have a solution of 
rank i+1 with respect to x. Then the chain must terminate after 
i operations, so that the system (E,) will be integrable. First 
we choose arbitrarily matrices A, and B 1 non-singular, of order n. 
Then the equation. 


(50) H = Ay + Bay - Cy = (0) 


will determine C,. Then a partieular integral is given by (43). 


The value of the K,; invariant for (E,) is determined by 


(51) Ky =~ Tae, Pe 
and then the relations (17) and (21) will permit us to calculate 


the invariants for the systems (E;.}),-.-,(E), 


Vize Hs 21 = Key 
By.) = By, 
-1 . wh 
Agey = Hyiz (Ar Cb Bi-vy Baea Mana: 
2 = oa 7. cl air 
(52) =H, A, yy + By Himlyiy 
Key = Ky, - Be Ag + AgaBycy SPaeata-2 


= ely _ -1_ “1 + 
+ B,_,(K,A,_jK, 77 il Ky) - (KyAs Ky iy Ky) By 14 


1: > =] 
+ KAR > BK De, 
These relations also permit the calculation of the coefficients 
A, B, and C. From the expression in (52) for Ay and B,_1, we 


can show that 





*, =1 , 
A = (Hy Bao »#eHjH) ‘As(Hs Hy pee: HH) + 


ce 1 - 1 
+ sa 
a (TT 4 (aah (TT. nie Ys 
B= a 


and, having determined K, 


C= By t AB -XK, 


Having thus determined all of the ifvariant matrices H Of); <i, 


“ly 
We may then compute the general integral for U by the iteration 
process of paragraph D. 

We may determine in a similar manner all systems which ter- 


minate in both senaes, and hence admit a particular solution of the 


form 


Be Weta: & co. er4 202 Gy x) 

which contains no integral sign. The preceding expression has rank 
itl with respect to x and rank jtl with respect to y. The sum 
i+ j is called the characteristic number of thse equation. This 
number does not change wpon successive applications of the Laplace 
Method. In passing from system (E) to system (Ey ) for example, 
the number i is diminished by h, but the number j is increased 
by the same amount; the sum of the two is unchanged. This is apparent 
if we consider the system (E,), with invariant H,®0. This is of 


the form 
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m 2° ( ooh tay U,)+B, (Oh + ag) = (0), 


which admits a solution of rank 1 with respect to x, and of rank 
itj+t+l with respect to y. For brevity, let n=i¢+ Jj and con- 


sider the substitution 


=] 
: ([A jay) 
where co is a colum matrix of unknown functions. 


Then gu a ss 
_ (Tsay) ™ — (Tasay) "© , 


Sy pAyd, = (ajay) 19 
Wy ay 


Using this (53) becomes 


(54) é 4 (Taar™ 2S} +3, { (T ayay)™ * = (0) , 


An integrating factor for (54) is miltiplication from the left by 


| Biax. Thus we may write 
+ JQ) = 
(9 fT aax (Tas 280 = (0), 


Denoting [2s ‘ (Jaay)™ = at (55) becomes 





which may be immediately integrated to give 


= . 
where Yy » designates a colwm matrix of arbitrary functions of 


y, and X a colum matrix of arbitrary functions of x. We know, 


already,that (56) admits a solution of the forn, 
(n) 
57) () = xX, + ®R yx t Bir + 20 Bn¥ 3 


in which @, @,..., @, are certain well-defined matrices of functions 
of x and y, and Y is a colum matrix of functions of y only. 

If we substitute (57) in (56) and take an arbitrary numerical 
value x, We find that Y and Yl are related by an expression of 


the form 


( (n+ 1) 
= + 
(58) YT Ay + At + e900 \ ogy ¥ . 


in which A , A por N are certain square matrices of func- 


n+l 
tions of y only. Substituting the expressions for @ and Y, given 


py (57) ana (58) into the first equation of (56) implies that 


(99) 2 {6 +8,¥'4 athe F 6.x" = “A(AYt A,z't. tA ae 


In order that (59) be a valid equation, the coefficients of like order 
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derivatives on each side must be equal. In this manner we obtain 


the system 


(60) I8-= 
0 ge “KA 


rd 


g8n + 6, = ANy 
@, = AN\nay | 


If we eliminate the 6, from (60), we obtain the equation 
hel 
(1) A ~ 2 (ory) + 2 (A Ng )— + +E) a a A\n)= Os 


a linear equation of order n+l with respect to Q{ .,.whose matrix 
coefficients contain functions of y only. If we solve (60) for 


the @ 3? we obtain the values 


(62) ©, atk Mn ’ 
-B, _=HKA, = 2M ni2) 5 
B= cK - 9 A T caret 4)j2 On oc N\ 
AL” 2, (oe nite 28 (oe aay, 
which permits the determination of the expression for @ . 
The relation (58) , between the matrices Y and Y, » per- 


mits an arbitrary choice of either one or the other for the defini- 


tion of the value of ,» by either (57) or the second equation 
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of (56). The choice of the former, however, will give the most general 
farm of the integral for aq , 

As an application of the preceding discussion, which is com- 
plementary to our original problem, we see that it is sufficient to 
choose n+e square matricss of arbitrary functions of y, namely 
N, \ eee oe] » S80lve the resulting linear equation of 
n+1°" order in one independent variable, (61), for the matrix CC , 
in order to determine the form of the equation (E;). The formulae 
(62) and (57) then provide the general integral of the equation (56) 
and hence (53). Repeated application of the substitutions of Laplace 
will enable us to determine the system (E) which has this general 
integral. 


st 
Now consider an arbitrary linear equation of (n+1) order of 


the form 
/ (n+1) 
u n+l ? 
where M, M- oo Mal are any n+2 square matrices of known functions 


y and OC .is.a square matrix ef funotions of y . Let a; Ne ni} 


be defined by the identity 


7 to i (n+1) _ 
a ntl +1 
eo) A =< MW = 2 (M,W) + oe (MW) “seo t(-l) a (MW) 


where W is any arbitrary square matrix of functions of y. (6k) 


then determines the /\'s in terms of the M's, for example 





hg 


N\= M - 0 Mh + aM "= #000 + _ jotem ae 


qo EE 


o 4 


(65) (23) el 


a1 


A 
nal 


In light of equations (60) and (65) we my set 


>. n+ 1 =L 
Moe. = 6-2) “AK. &, 


That is, if OC is a solution of (63), then we may determine the 
functions A j by the identity» (64), and the G, by the system 
(60). Upon substitution into (57), we see that the general integral 
contains the colum matrix Y and its derivatives up to order n,. 
Thus we have proved the following: 


Theorem IV: In order that the linear equation 


ao 2¢ | ate) 


atl with respect to y, that is for the 





admit a solution of rank _ 
general solution to be of the form 


@©= x+@ ¥+@,¥4 on +B yi) 
4 





such that the solution cannot be written in analogous form in which there 


appear fewer derivatives af the arbitrary columm matrix y , it is 









necessary and sufficient that AH, 









a linear equation of order nm +i. whose coefficients are functions of y, 


that G does not satisfy 2 simtlar linear equation of lesser order. 





Jo The results given in the section have been kept in close analogy to 


those contained in Chapter II, volume 2 of Lecons Sur Le Theorie Generals 
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Des Surfaces by Gaston Dsrboux. In this chapter Derboux accomplishes 
more With the single second order linear hyperbolic equation than is 
available to us with a system of such equations, due to such causes 

as the non-commitativity of the ring of nxn matrices over the field 
of functions of two real variables. It is of course probable that some 
of the results obtained by Darboux which were not considered by us in 
this section, could be obtained for systems by matrix methods. The 
reader is referred to this very fine work by Darboux which has been the 


inspiration for much of this theaig. 


+ 





ory 
mult ‘eed we a. 


A. Let us now turn onr attomtion te + 


g 


order in thrsas indepemdent varishl 


(1) ne (uj = es xs ee ~ ou cr 
tay as Pd ‘- ae 
MmeecoellicientS 2, &, c, Gd, & FF, and 
s 
Functions of x, ¥, and 2, combinnouels 


wo my need. Wa wish ts stouaecsxs Gals 


‘4 ~y 


arabe 


eu 
Se ath 


cussed previously, 
Pirst order equation 


the equations, in the hone tbat ato. 


chain will terminate with vwanlarire 


’ a = ok 
macaged AB Paez corie 


system can thon be 


To commence this opecauti. 


different types, naroly 


(2) uJ = uw + 
(3) we = Yu ios + 
(4) i ee 


Let us stervt with tne change of veriaslea 


the operator h. whies appasre in (1), 


INDELENDAM? 


LO LY 


CF HIRD ORDER IN 


Toes lew i = § 
L£< a 2. CL 


«ug will then attempt to cascade 


“4 
rina. ad deme we) 


bu 


O 


WH POE 


a 


iL) 


reaucing { 


). 


8 
aS. 


VARTABLES 


a 
are 


pare: 
NY fabs: 


on 
I, 


he linear ogqnations of third 
as whien have the form 
24 + Gly 4 Pa, + gu 


ae to be considered as 


.6© 28 many times as 


- 
—_ 


O, 


1 the manner of Laplace dis- 


‘8 
4, 


to a system of three 


cu 
J 


Cesived . 


iwaper of iterationa the 


Suvariantu, and that the resulting 


+ can, 


im, We wat consider exbstitutions of thres 


Usine the definition of 


Lae 
¥ 


verify that 
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(5) a + by + om + au K(u) + kw, + Lu, + mn 


where k = a + ae - 6 


al Z. 3 
= a + ab o- f 
(6) lle ? 
= a + ea + dba + 2a = 
“al Jz y Z yZ - 


are defined as the first three x-invariants. It is apparent from ths 


symmetry of the equation that we could have just as well made the substi- 


tutions ; 
ie) wy = wy + bu or 
(2° U3 > wy + cu. 


These lead to equations similar to (5) with the y-invariants (arising 
from (2')) being 


hy = > mts be . d , 
Lo = » + ba -~- f, 
n, = ryt cb, + ab. + be = Ves 


and the z-invariants (arising from (2'')) being 


0 =x ¢ eb - 4a 
3 eo 3 

Kk, = Re + ca = e, 

m <= ec + be + ae + cf = & 
3 hy x ay e 


It is also apparent that whatever results we obtain from substitution (2) 





a5 


would be symmetrically obtained from substitutions (@') and (2''), 
hence at this point we will limit the majority of the discussion to the 
effects of substitution (2). 

Should we be rewarded by finding k,= bh, =m = 0, then (5) 


will become a second order hyperbolie equation in only two variables 


(7 + buy + cu + due=™ oO, 
) Palle ii 1, =e 


We may then apply the Laplace method to (7), as described briefly in 


18 
Section I, or in more detail in the works of Darboux. 


Consider next the change of variables (3), which leads to the 
equation 
2 2 2 
(8) uy, + euj = K (u) + how + Kyu, + |e, + mu, 
where now ms = a, + (ab), + Cay + abc -«- g. As before, the 
symmetry of the equation indicates that the following changesof variables 


would lead to equations similar to (8): 


v4, + th 


us = % + au, , 

ty =u, + cu, 
(3') > 

1 

u, = * + au, 


These substitutions lead to the respective invariants 





5 


He 
tt 
a? 
N 


+ (ac) Te ba + abe - g 


of Z ) 
m~ => Of tay = Fe 
2” “xz Zz —. Ps 
35 mm 
m: = »b + (b + ab ‘b - 
- (be) } + abc B, 


m = ec + oo + be + abe - g 


es ce + (be) + ac + abe - g 
3 xy iG J 


1 = = 2 == x. 
In the case of equation (8), if Kk, Ll h, m 0, we will 


have reduced our system immediately to three first order equations 


Eo - an = Uy ‘ 
+ = # 
(9) Oy M1 El, 


J 
2 


os 
i. 4 eu 0 
Which can be solved in inverse order by quadratures. 


Thirdly we consider the change of variables (4), which leads to 


the equation 
(10) er + au, = a (uw) + Lou, + ku, +My wu. 


where mj,= d, + ad-g is the first -x-invariant. Once again 


symmetry consideratiams show that the change of variables 


urn = Uys + avy, + cu, + eu , 
and 


43 = yy + aan, 4+ bu, + fu, 
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lead to the ey-invariant and -z-invariant 


Mo ey + eb - 8, 
m3 = ry + ime - &. 
respectively. 
If we are fortunate once more, such that i, — x, - == 0; 


we may solve (10) by quadratures, and then apply the Laplace method 


to the non-homogeneous equation (4). 


B. We have in the discussion of paragraph A encountered some eighteen 
expressions which we have labeled “invariants". Organizing these in 


somewhat more orderly fashion, we have the following table: 


+ be - -d 


ne 
it 
Ps 


ha = ec. + be - @ 


], = > + ab - f 


m = a + ba, 4 Ce) eo Se 


+ + = 
tT ab, cb. be g& 


nv 
tl 
io 

Ps 

N 


m ™ Cc + ac + en HCE - 
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my, = ad t ad - g& 


lf 


2 sr he = @ 


ms = £4 > 67 - Bg 


me = &, + (ab), + cay, +t abe = g 
3 = a 
ot Sys a rl), + ba, + abe g 
1 
= + = 
mn, a 4+ (ab) eb + abe g 
m = v. 4 (be) + ad + abe - g 


3 


e +(ac) +be + abe - g 
xy o tak 


il 


wo ues 


e +(be) tac t abe - g 
xy Ray 


We now wish to investigate the character of these eighteen “invariants”, 
as was done in paragraph B of Section III for systems, in order to see 
if the term “invariant” is appropriate. That is, we are concerned with 


the change in these sxpressions when we make the various changes of 


variables 

(11) uv = ZX u, where rN = \ (x,y, 2) is at least 
twice continuously differentiable} 
t 

(12) Se as — (x dg 


I = W (y') 
z= Liz); 





at 


Let us consider first the change of variables (11). Equation 


(1) is then transformed into a new equation of the same type 
yt t 
(14) XK(u) = uy, teu, + dy, tewy tau + eust 


+ fu, + g'u’ =0 


but with new coefficients 


g lies 
a eae 
bi = bt d SA 
to 1 Or 
oo 2 oh be : 
i Ws. bo, Et Cw 
(15) SS ON Uy noe X yee 
os te | ty) © Ang io a a= 
vA dx oe >» x92 
f= fed Dyed ,1 QA 
» oz tor * RK OBS 
| 2 a4 Daa ore a°A ay 02> os 
Sie et He Sete gh ote A Seay tage 


Using these new coefficients, let us compute the three new x-invariant, 


k,, |, and m. We find 


Cre) ky - a. rac «= + & Ki, 

% al é tte ' = 
(16) 1, a + ad fe ,, 

t t Le | e ¢ tt 

m= Ant te + oe + ie 9 ae ele 1a od , 
The new - x-invariant, m4, becomes 
; t t '_t ' k » | or 
(17) ma tad - ¢ fe 2 > Ge d 1 me. 

re . or Oe 


g 
While the new me invariant, ms ,» becomes 
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ms = ay + (a'b’), + c'a' + a'd’c'- g' * 


oy 
_ = ma % Dd SAHA, 
hae YI NSE 

From gyumetry, similsr resuits are obtained for the other 
h', k, 1', andm invarients. Thus we wee that the k, h, and | 
invariants, which appear as coefficients ef the x, y, and z derivatives 
of wu, respectively, are indeed tive inveriants under the change of 
variables (11). The m invariants, on the other hand, do not repro- 
duce themselves exactly under this chenge of variables, but instead re- 
produce themselves plus & linear corbimation of the true invariants. 
Hence we shall in the future refer to theses m invariants as quasi- 
invariants. (These will appear again in Section V.) 

If we should make the changes cf variables (12) and compute the 


resulting invariants, wo will Yind thet 


ky , D Zot Ky 
PV], 
(19) at = Ql. 2: m, 


, ? 
nm Vos Yoyo Z gs ms 


3 


bh 


3 


i 


mye Ox Ppl 


Finally we note that (13) merely ees the x-invariants 
with the y-inveriants, the y-invariants with the z-invarients, and the 


zgeinvariants with the x-invariants. 


C. Let us now develop the cascade of equations in the manner of Laplace 


and Darboux. Consider first equation (5) im the event that at least ons df 
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the invariants is not zero. We wish to transform (5) into an equation 
for u, which is in the form (1) , (See in this conmection remark in 
footnote (5)) In order to do so, we proceed to integrate equation (2), 


obtaining - 


(20) s 


Saas | . (eax 


where X(y,z) is an arbitrary function of y and z only. Then differ- 


entiating (20) with respect to y and z yields 


ie 2% —-_ fol ayes +3 - 


J 
Sb Tes 


—.., om 2 Sy ee ; _ 2. (\ (eax) es Sect 


Substituting (20), (21), and (22) into (5) we obtain 


(21) 


u, + bt, + ae + dy = Me ~ {pax ee J see he: + :|* 
(23) eae (eax 2 {as J x ; 2 ( § eax) (a ; cl 
me, -aax 2 i ‘a , rt - a (Ss) f (0 ujdx + “} 


& 


Now (23) is an equation for u, alone, but it is an integral- 


differential equation, interlaced with arbitrary functions and their partial 
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derivatives, and hardly in a manageable state. Our aim is to reduce this 
to an equation for u, in the form of (1) and to this end we will dif- 
ferentiate (23) with respect to x. Before we do this, however, it will 
be necessary to make certain assumptions regarding the invariants and the 
coefficients. If these assumptions are made, the integral terms and the 
arbitrary functions will disappear when we differentiate. We will con- 


sider a number of different methods. 
METHOD 1 
The first assumptions in this method are 
(s) x= 1, = mF O; 
(b) a isa function of (x) aly. 


Under these hypotheses, (23) takes the form u a bu, + du, = 


Ly 


4 _ less (09° ayes +x 4 2 [ona + x{ tof (a 4 x 


adx 
Multiply both sides by © » and differentiate with respect to x to 


ra 





get 


= r ( + bu, + cu, du) 
oe u u . a 


eh 
( >= _ Geax 


Ox hy 
( adx 


a 
m1 


(“leg + Lon + Shi. + du, ) + es 


+ ae + buy, + CUpgy + duty + ey + byUy ¢ 4,u,) = 


im - u, + ie u, + e fea, 


L 


Z 
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- \adx 
Now we may multiply both sides of (24) by mye S and collect 


terms to obtain 


+(a- Qi 
_ eae gy a eels + an + sail + 


+ (c(a - dlog my ) + C, - Mm, ) me + bla - dlog My ) + by - my ) Uy, + 
Z 


x x 


+ d(a - dlog m) 4 4, - m, ) Uy = oO, 
x 


Designate the coefficients as follows: 


a; ™ a- dlog m, , by = b 
ox 


ee oT 7 Gan - Mm 5 f, = bya) + > - i 5 


cj; =c d,=d 


) ? ) 


Then (25) is in the same formas (1): 
Our new x-invariants are, fram the substitution YW, = ie + 4,U,, 
ky, = ayo + “ = Oy = 4 7Cy + As - ayer - $3, 7 a = 


= 41, = 2, ¥ aL< 


4 
my = log My = Cy. 
yioz 


Thus 


k, =k, ~ (cy + Frog m \ 
Szox?. 





Also 


2 a é 
I= 212 + 1, - f£, = 4&)b)- 0 *10¢ Mm, ~ &)b) +, - i= 


oyox 
L-(», + )*108 m ) ; 
0 20x 
my, = a,d, + a1, + “1"Ly + “99 = 6) = 


ee” 5 - = eee d* log m, 2 ee 910g m, 
020% OKO 


= 9° tog m 
OxdyO2 


= 2 2 2 3 
= m f= b oO log m, 4¢ a log Mm 4 Q log 2, 


; ya Jyoz oz0732° 


If now the new x-invariants are all zero, we may proceed to reduce 
the equation as described in paragraph A and solve. If, however, at 
least one of the invariants is not zero, we are in a dilemma, for the 
hypothesis with which we transformed (5) into (25) are no longer valid in 
the general case. We can no longer continue our chain. To enable us to do 


50, We mist change our original bypothesis to the following: 


(2) x, = l=, # 0, 


(b) All coefficients are functiongof x only > 


tt 
tt 


(c) b d 


= @ 


Under these hypotheses, our invariants are also functions of x aly. We 
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may transform (5) into (25) as before, and the coefficients will have 


the same form. The x-invariants will be different, however, and in fact 


Will be 
ky = ky a Cy = ky = oS 6 
(26) Lyon), We, = OP. S, 
ay x 


Hence k, = |), = m,, our hypothesis are again satisfied, and the chain 
may be continued. 
Observe that, as our chain progresses, our x-invariants will always 


have the form 


k ja. 


- JSC, Oe n:e 


14 3g 7a4337 agg? 


Thus, if our invariants are not originally all:zerg gthey can only become 
zero if they are an integer mitiple of the partial derivative of d with 
respect to x. 

Note also, that to have a chain of y-invariants, or z-invariants, 


our hypotheses (2), (b), and (c) must be modified accordingly. 
METHOD 2 
In this method we assume first that 
(2) k= 1, = 0, 


Under this hypothesis, equation (23) takes the form 
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“20x adx 
(27) Mog. t+ Pe) + Mg + duy = — is uyex + | 


adx 
Multiply both sides by oa and differentiate with respect to 


ser" 
x to get 
(adx 
a (4. - d log m, )Nalegee bu, + cu, + du,) + 
™ a f= 
adx 


1. _ + Pa. + rs e oe pli wi. ae ce 4 du, ) = 


my + 


-\adx 
Then multiplying by m) a 5 and collecting terms, we obtain 


2 + (a - 3 Log m3 ) Ble ait PF + ones +. ae. + 


gx 
(28) + l° (a - dO log my ‘| “ie + |b (a - d log my )+¥s | Uy, + 
“Ox. ™S : 


a SAX) se - ee, 


As before, designate the coefficients of (28) as follows: 


= ae Od log m, ; by = db, Cy = Cc 4 a,™ d, 


= 


eal 


Then (28) is in the samo formas (1), and is identical with (25). 
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Now, however, the x-invariants resulting from the substitution 
uw, = Yh, + &)U, are 
ky = aye, + ay 67 Saye, 87e) + Se ae 


=a - ¢ o 4102 my 


dz0x 


é m. 
] = a é = Sn oo O tog 1 
h 7b) “ly fy ay x es 5 
™=aid, + bia ca ~ = 
Bie Fe2ity + 01% + Seay + Goel sy 
= + ba + ee a 
ma - dy, a ca, + Sy, ~ og mm, .. 
Ozox 
3 
~c 0*log m .o logm 
a 99x dOYDx * 
if kh = 1), Sm), = 0, we may reduce as before, and solve. 
If, however, at least one of these invariants is not zero, we are again 
faced with the dilemma that our original hypothesis is not satisfied in 
the most general case, and the chain cannot be continued. Thus, we mst 
again modify the hypothesis in order to iterate and continue the chain. 


Our hypotheses becane 


(b) a =se(x), = b(Yyz) , o = C(Y 2) 


a= d(x), oe = e(x,y,z)=ac, ff =f(x,y,z) = ab 4 


g = g(x). 
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Under these hypotheses, =ad - g, is a function of x only, as is 
O log m, Now when we transform (1) into (25), the coefficients 
QO x 
will be 


faa eo - dlog m 9 ny = 2b, o='& 5 dj = a, 


> 


op 275)» f) = bpe,, 8) = 048, dy, - me 


Thus we will find thet 


4 


k, = 0, yy =O, m =m-gd, 
hence, the hypotheses (a) and (b) are both satisfied once more, and 
the chain may be continued. 
Observe that, as before, as the chain progresses, the invariants 


must have the form 
Pa By = te, 0 3 3,0 a ~id., 3 i= Oe 12 acs o 


Thus if m, ¥ 0, the invariants for the je iteration will vanish only 
if m),- is an integer mitiple of d,, 1.e., my ajcd,, for same j=0,1,2... 


We also remark again, that..to have a chain of y-invariants or z-invariants, 
hypotheses (a) and (b) must be modified accordingly. 


D. Now we wish to consider methods for cascading the equations when we 
employ the substitution (4), in the event that at least one of the 
ex-invariants is not identically zero. For this chain we wish to solve 
(10) for u interms of u_j, and, after the prescribed differentiatias, 
substitute into (4) to obtain an equation of the form (25) in u_j. 


In order to accomplish this we are again forced to place certain conditions 
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on the coefficients and the related invariants. Thus we introduce 


Method 3: For this we initially assume that 


(a) l, ~ k = 0 e 
Solving (10) for u, we obtain 


i 
mae) 


Lu= 


Cea + Gu. 7 ) 


Uy = 1 Maat ae + ayu_,)- zo — 
=» 1 (u, + au +a,85) -1 _, @_, +08 _. ) 
we 1, © 3,83 2 oak je Med’ 


i 
4 
B 


. #2 OQ m3 (s ® x5 
= ee * lp €4,u ,)-k gk (u_, ‘tau. + ayu ) 
~1 d 





2 
i. omy 2 om, m 
a. 2 5 Dz (uly + wr en — 0 “= (uy ae au_j). 


Substituting these expressions into (4) and collecting terms we obtain 


(29) Sijcs 


b 
2, i P21 Liye + 1 =1y, + © 1 -Ley - CW at Pot, + 
+ le “ 8.144 is Or, 


where &.,= 4, b4= db - 0 log My » Cy Ce - do mat ; 
ay or 





68 


"7 2 a-@ d log mj: b 0 log mi: 10g mo1 


= 0 2 Ovo2 


+ d log Mj 0 log m_} 
dy oe) : 


$.)™ Pailic t S221 3 fa ~ eal, +a_jbiy 3 and 
S.1= Gel, + b.ia_4 “19-1, +3505 - moi 3 


and (29) is of the desired form. Computing the new -x-invariants, we 
obtain 
2 
= - - 3 log m_ 
15 ay 2 ety 21s 
Ox9d 
k =|€e = a ~ Ope di. 
6 x 2 ~ 4 
a xO2 


and 


fo mi4+ dy - a (c log mi + b d log Moy , 010g oe 


= “ot Ogos 


d log m7 O log es.) me, ay (c = 0 log ms) . 


oF de oF 


ye (b - Dice m1) + Ayo , 
J 
As in paragraph C, we will be in 2 dilemma if the new - x- 
invariants are not all zero, fkethen our initial hypothesis (a) will not 
be satisfied in general, by (29), and our chain cannot continue. To 
escape this trap, we once more impose additional restrictions on our coef- 


ficients. Our revised hypotheses will be 
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(b) a =al(x), »b =b(y,z)> =c(y,z) ; 


a =d(x), o =so(x,y,z) = a(x)ec(y,z) , 


f =f(x,y,z) = alx)bly,z) , « =elx). 


Under hypotheses (a) and (b), m_, becomes a function of x 


only, and hence the coefficients of (29) become 


ai4™ 4, dy = Dd, c.. = €. doy =, Out = 6, oy = 


With these coefficients, the -x-invariants for (29) are 
PS | 3 _ x 
1. = oO Oren 1m ee: 
which means our hypotheses (2) and (b) are again satisfied, and the 


chain may be continued. As observed in methods 1, and 2, the invar- 


fant mj _ 3; can vanish only if m1 is an integer mltiple of d,, 
that is, m_}™ -jdy. Wo may add at this point, that i? m is a positive 


miltiple of d,, we should use substitution (2), while if it is a nega- 
tive multiple of d,, substitution (4) would be the more advantageous. 
One more it is worthwhile to note that if we desire a chain of - jy - 


invariants, or -z- inveriants, hypotheses (a) and (b) must be modified 


accordingly . 


E. Finally we consider a method for cascading the equations when we employ 
the substitution (3), in the event that at least one of the invariants 


in (8) is not zero. Although this substitution appears the most natural 





(0 


in that our equation (1) is immediately broken down into three first 
order equations, we will find thet the most stringent conditions on the 


coefficients are required in this case in order to generate the chain. 


Method 4: Our initial assumption in this method is that 


(a) ho=k, =l, = 0. We must consider the system of equations 


wy + au Uy 
—_s 


= + he BH» 


2 ——e 
ie + “= mu , 


3 


If we solve the second equation for uj, substitute this value 
into the first equation and solve for u, and then substitute the result- 


ing expression into the third equation, we obtain 


2 2 2 ies joan [ay 4 


u + cui = m 
i 1" le l 


+Y (2) f dx + x0)| 5 


where Y and X are arbitrary functions of their respective arguments. 
We then proceed as before: 


Se fea] apn igen a 


1, 
a 
z 





Taking the partial derivative of both sides with respect to x, 


(adx 
~) te - dlog my (ur, + cu, yeus teus de “| 
an = yall lyz 1. wl 
—— \pay .: 

oh fe 


uj»dy + ¥ 





GL 


(pay -aax 
We then multiply both sides of (30) by «6 and take the partial 


derivative with respect to y, to obtain 


| (»- OLogn=) (= ~ DroaZ (uf teuf)+ uy tem, + eit + 


= 2 : 2 2 ; 2,, 2 2 2 
+(a, - @ Logm) )(u;_ + cuy) 4 (a - dlogm;) $a oh ony + Cyt) + 
0 <XgY 2 
tS gat Magy + °y%q + eal, + “| _ 2 bay o 
2 \-bay 


Finally, multiplying (31) by ue » and collecting terms, 
we obtain an equation far ug, which is of the same formas (1), but 


2 i 2 . a 
whose coefficients are a, =& - Q Leen, ; by 2b = OLogmy ; ce] = ¢; 
Wer oy 


ee), 2 SS 22 ie eC, 
ceed ys) 681 £0) En i = appt =~ ? 


eHsey 1reqbyes +. (aye)y Hey ~ my s 


Computing the new invariants, we Pind 


_ §2 2% 2 ee 2 
he = by, + biel - gj = >), -¢1,, 
2 aD a 
ie «8p, at Ga") -) a, 7S 
weap ter tf Eo, 
m = a4 + (agb7), + eat, + acbics - aS 


2 2 2 onze 2 
= my, + Lee + (azb3), a at Cy 


AY © 


We are faced with the same difficulty as before. In general, the 
invariants dass and ky, will not vanish, hence it will be impossible to 
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continue the chain in this manner. We will be able to continue the chain 
if we revise our hypothesis to be (a) ho = ky =1,=0 5 
(b) All coefficients are functions of 
y only, except c, and that c 
is a constant, 


With these hypotheses our new coefficients became 


2 2 —_ e > 
a, = 4; by b ~ d Logm 5 cae és a° = pies 3 
2 2 2 . 21,2 

° OLE ee Sere ae 

6, = acy; fy aybr+ hes gi =a,bicey - m » and the 


new invariants became 


m= m + ca 


Thus the hypotheses (a) and (b) are satisfied once more, and the chain 
may be continued. We observe that in this mothod the chain will terminate 
with all tnvariants zero if and only if ms is a negative integer multiple 
of Cay. 

These methods described in paragraphs C, D, and E are not,of 
course, the only methods by which chains of equations could be generated. 


We could do a method similar to method 1, say, with k, =O and a Fe, 


a 
or we could do a method similar to method e with ky the non-zero invariant, 
or ly the non-zero invariant. There are many variations on the theme 
presented here, but it is the theme itself which is the important thing. 

In any such method certain restrictions of the coefficients are necessary 


to iterate once, and other restrictions are necessary to iterate more than 
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once. The basic idea is to find an expression for u in terms of the new 
dependent variable, which may or may not be an integral expression. We 
then differentiate this expression as required by the particular equation, 
L.e. the equation involving the invariant coefficients. Substituting the 
resulting expressions into this equation, we reduce it to an equation in 
the new dependamt variable alone; if necessary, we separate these ex- 
pressions fran the invariant coefficients and differentiate again suf- 
ficiently to remove any integral signs which may appear. If we can 


accamplish this, we can in general iterate the equation. 


F. We consider now the most general form of (1) which can be reduced by 
these methods described in paragraphs C, D, and E. The question is, what 
mist the coefficients of (1) be, in order that (1) may be reduced to a 
second order equation by these methods, when the x-invariants or the «x- 
invariants are not all originally zero? We observe first, in methods l, 2, 
and 3, that in order for the m-invariants to vanish after n-1l iterations 
we must have 


g - ad = nad, , where n=l, 22, ... 


(32) d+ 


d= 
n 


|e 


Equation (32) is integrable, and integration gives us 
Sa { 
(33) ds e e g dx plus an arbitrary constant 
n 


which we choose to be zero. 


Using (33) and the hypotheses of method 1, we see that our original 


equation mist be 





(34) 
4 e(x)(H+¥ +42) 0, 22,2, core 


for method 1 to be applicable,’ Using (33) and the hypotheses of methods 


& 


2 and 3, we see that ow original equation mist be 


Sara e a(x)uy, + bly,z)w, + e(y,2)U,, +( = ax t 


(Bae ¢ \Rax 
\- 


(35) 
+ a(x)e(y,z)U, ta(x)b(y,z)U, + e(x)Y = 0, RB =1,2,000 


3 


for methods 2 or 3 to be applicable. 
Finally we consider the form of the coefficients when method 4 will 
be applicable. In this case we observe that in order for the m invariant 


to vanish after n-1 iterations, we mist have 


g 9 abe <a. 9 n= Lye, 000 
(36) ap, ba=g, 
nm 


Equation (36) may be integrated to give 
(0 
(37) a= 6 a \s* gay . 
en 


Using (37) and the hypotheses of method 4, we see that our original equatim 


must be 


AY 
~ (Bey (27 
+ eb(y)u, + (ce n ‘. nm aa)\u,t 
cn « 


flee 


- (bay {bey oe 
wate 5 {° a -aSy)u, + (yu teu, + 


(38) 





+ (2, + EDe —_ (+ \2 


for method 4 to be applicable. 


2 


Cc 


n 


ey } u, + e(y)u ao 


{E 





SECTION ¥ 


HYPERBOLIC EQUATIONS oF N™= oRDER 
IN N INDEPENDENT VARIABLES 


A. Having considered the extension of the Laplace cascade method to the 
equation of third order in three veriables of the mixed derivatives type, 
we will now turn our attention to a generalization of these results. That 
is, we consider an equation of nth order in n independent variables, 
in which only the mixed derivatives appear, and we ahall attempt to show 
how these methods can be applied to such an equation. 

Before proceeding to the ne2 order equation we will first intro- 
duce an operator notation which will simplify samewhat the invariant ex- 
pressions and relations, which tend to become awkward as we increase n. 
in Seier to do.sq, we mention here a rule which has been devised by J. Le 
Roux’? and shall hereafter be referred to as Le Roux's rule; 

"For any differential operator D which is a polynomial in Q 


sx 


the following relation holds; 


Cre 


ena 0 
oy 


x aye 


D(uv) =u0(v) + gu D,.(v) 1. ou Di(w)t Bs a Diy (v)t2 a Day (¥) +. den »" + 


ar 








Or 


d>u p™ () + 3 33a p(w) +3 23u Dywy(¥) + d3u 9" foee 


Q x3 a= dx" OY axdy- ay3 


¢ 


where D dp » Dz dD ’ ete. 


onus 


ie) 4) 
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ie 


This rule is an extension of Taylor's formula, and is, of course, 
canpletely analogous to Leibnitz' Rule for the nvh partial derivative of 
a product uv. It can be similarly extended to produets of more than two 
functions. As an example of this rule, consider the operator D(u) = 

su_tau.+bu.+cu =0. Here D = gd + ad 4 bd + ¢ , hence 

X 
xy v dz dy ° Ox Oy 


g 
bd, = 0 ca , Dy = 0 + Dio Dey = I, the identity operator. 
oy dx 
For this operator D, it is a simple matter to verify the identities 


(1) DD (u) = D(u) + hu , 


(2) DDi(u) = D(u) + ku 


) 
where h=a,tab-c, k* by 4+ ab-c are the two Darboux invariants. 


Then denoting vu, =D;,(u) and observing that D(u) = 0, we have frem (1) 


(3) Dy(u,) = hu 


e 
Operating on (3) with BE, we obtain 
' ¢ t 
(KY Dipt(u,) = Dinu) , 
Using the identity (2) to evaluate the left-hand side of (4) and Le Reur's 


rule to evaluate the right-hand side, we find 


D(u,) + ku, =h Di(u) + u Dit (u) 4 ot Dy (u) | ere 


hu, + Jhu 
J 


hu 4 dhe bh wu 
1 ae A 


J 
h Ww, + ee D(a), 


il 


or 
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D(u,) - do D. (uj) + (k - h)u, = 0, 
y 


This may be written in the form D, (a4) = 0 


where Dy = Q 0 «} a7 Q =p by 0 + CL iy 
yz oy 


ox oy 
ay 
a > 


ec, =c - b dlo +b, - a; 
y 
and these coefficients agree exactly with those obtained by Darboux”° 
for the first equation of the cascade method. 
Our problem now is to expand these procedures to the n*h order 
equation in n independent variables, of the mixed derivative type. Can 
we in this manner predict the form of the "invariants", the corresponding 


identities, and methods for cascading the equations? 


B. Let us, for convenience sake, change the notation slightly. Instead 

of indicating the independent variables by x, y, Z,ecco, Wo Will find it 
th 

more advantageous to label them x), Xpy 0eeeX,0 Then let the n order 

operator with which we will be comeerned be defined by 

(5) Due Day 


u : i,=0, 1, 2 n 
i i le ih Maa 
‘n? ) 


pd, ito" 


where we always require that i, = Aouad and that 1;5 1,,] if and only if 
1441400 Tho sum is to be taken over all acceptable camb inat ions of the 


i, where i,50,.00,m0 In (5) wo write u,= 23 .» JO, ly coon: amd define 
J d J ony 
ou = Uo Finally we require that in all of these operators we will set 
a o-8 on? u 





ie 


Qin °°? yn =1. To illustrate, let n=2. Then 


D(u) = 2 ey + &) ou + Ago Ou + Ampile 


For n= 3 We heave 
D(u) 29 3y ten, Ie ye, FU +8 OU 
apasty °° gaps,” Sor gs, “or ge,* 


+ So) OM + So02 9% + 8903 IB +2qqq" » 
O* d¥, 2 


% 


and the brevity of (5) becomes more apparent with increasing n. Having 
geen these examples in explanation, there will be no loss of clarity if 
we drop the zero subscripts of the a's wherever they my appear. Thus 
for n =2, we write D(u) = Wo + Spy + ay + aus 

In order to point the way, let us briefly summarize the end order 
results of Darboux and ow gra order results in the light of this notation, 


¥oxr the second order equation we have the following operators: 


D=e QQ O_ a, 9 a, Q a 
a, * 7 + Ve ee + 


dt gXo ? 

ft 
D= a 3 

L = 4 

d*, 

D. 

= +a 

tf 
ae 


The two Darboux invariants are given by 
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t 
= —— D = 
a, + 2,8,- 4 28) a 


0 
oxy 

k= 9 8% $a,4,> 9 = Di(a) - 4, 
3% 


while the identities (1) and (2) becoms, respectively 
D3D, (u) = PD(iu)+ hu =D(u) + { Da(ay) = a u . 
D,Do(u) = Diu) + ku =D(u) + { D3 (a2) - af u. 


For the 3°¢ order equation, the operators are as follows: s 


D=d Q 0 +a d do 4a d oO Q + 
$x J7%2 0x3 ou BS pte gs 


a 2 + a 3 a O a 
ee ge 
D=d a a e a , a 
i jndx, + “fet ae le 
pled Q a a a, 
ete dey 20K t Bie * FD 


ee P 4 

12 jx, * 12 » 

mg e+ 3, 

Dag = 9 + 903, 
=) 


and associated with this equation are the invariants 
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"! 2 mY 
hy = Dy5(a)2) 2, mn = Dy oD, 3$854)-2 
a tf 3 e tf 
B= Big’ 12) 2 ee ogy 
1 
= Dip(@p3)-8, a oy Dy, (2, )-a 
_ me 
5 > Do3(8 5) -B, m5 = D5, Dip(a,4)-2 
tf 1 
1, = Dy 3 (253) -a, m3 = Di Digl 20 ) =a 
L, = 0, 3(8,,)-, my = DAP, (%),) 
t; - q 
m, = Dd, (a,.)-a ‘ n= D53(a,) -2 
m = Dz (a,3)-2 m= Dj x(8p)-3 
t 0 
mn. = D, (a,,)-a n= Dip (a3) a . 


A typical identity involving these invariants is 


DiD33(u) = D(a) + (Dyp(ap4)~ap)uy-+(D, (854) ~2,)¥, +(D, (a59)-a)u , 


C. Now let us generalize these results to the n’h order linear differ- 
ential equation of the form (5). In order to demonstrate the magnitude of 
the problem, let us first compute N(n), the number of different identities 
associated with the n’® order equation, D(u) = 0. 

An identity for D(u) is a cambination of differentiated operators 
p! ’ D ati. such that successive application of these primed operators to a 
function u = u(X), where E = (Xy%y000yXp)> will yleld D(u) plus a linear 
combination of the aired derivatives of u of order less than n. We shall 
denote such a linear canmbinatio by R, for remainder, and define the coef- 


ficients of R to be the invariant coefficients (or simply invariants) 
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associated with the operator D. We shall show presently that the order 
of R mist be less than nel, although this is not evident a priori. 


To be a systematic we proceed as follows: 
q 


Select arbitrarily a primed operator having order nel, i.e. D, 
for some i11,2,...,m. There are, of course, n choices for this oper- 
ator. To obtain an identity es defined above, there is only one choice 


f 
of operator p'?? such that D,D (u) +R. This operator must be 
(n-1) 
27 Ceoetel £41.s.n 
t 


i 


, since this contains the differentiation © which 
x. 


is missing from D Therefore we have exactly n identities of the farm 


t (n-1) ( aa : 
a Dd, eee tel i¢il. eer v) D(u) + R, 
L152, ccoyNe 


7 
Next select arbitrarily a primed operator of order n-2, D for 
aS 


Some 171,2,...n and some j2=1,2,...,n, jJ#i. There are n choices 


1 . 
=D. , we total only 
4 a 


m2 (ma) S 2) such operators, where (7) denotes the appropriate binomial 


for i, and nel choices for j, but since D 


coefficient. If we wish to obtain an identity when we apply such an oper- 


ator, wo must apply it to one of thresthings: 


(n-2) 
1) er 2 oee i=l] 141 eee jel jel te 


(n-1) p\a-2) (a) 


2 
Ds 2 woot Hl i+1..n L 2 eee jel jel eeell e 


(n-1) (n-1) 
3) Dy o eoogjrl jel Se Oo ecooktl itl em 


Therefore wo have 3() identities of the form 





Ds, L(u) = Diu) + R 


where L(u) is one of the three expression shown above. 


If we continue to count the identities in this manner, we find 


that 
nel 

(6) N(n)™ D> a, (3), 
gel 

where 


Mca) Case) Ge) 


Kt? the set of all ordered.integral partitions of Jj, 


(7) ‘> 
As 
of 


_ s 
cares J kell?, neh G- 
{sl | 


To verify this forma it is only necesary to count, as we have 


already done for the terms J=1 and j=2, the identities for the J” 
term. For each of these identities, the last operator to be applied is 


p\9) 4.9 where the i form all the possible combinations of n 
q 2 oo7°e 3 Dp 
integers taken j at atime. Hence there are ( 4) choice for this oper- 


ator. To complete an identity we mast apply this operator to a combina- 


tion of operators which contains each of the differentiations Q ; 
aA 
P=—1,2, oee, J, Once and only once. Thus, for example, we may have a 


combination of j operators, each of which has one differentiation 0 ; 


x4 
Pp 


or we may have a canmbination of j-1 operators, ome of which contains 
2 
the setond arder differentiation O 1, while the otherg ecmtain’ 
Pp 7d 
only first order differentiations 2 s and so forth. The number of 
a 
“?P 
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choices for each operator can only be the number of combinations of the 
differentiations still to be included, taken as many at a time as the 
order of the operator being considered. Finally, since the order in 
which these operators are displayed is a factor in determining which 
identity we have, we mist consider all ardered integral partitions of j. 
The summation in (6) is taken only up to n-l since for j=n, we have 
the trivial identity D(u)=D(u), and this is not to be combed in our 
determinations. If we tabulate the appropriate values for n=2,3,4, and 


5, we have 


s 


= 2)=e 
a, = 1 N(2) 
a N(3) = 12 
a= 413 (4) = Th 
3 


D. We turn now to a consideration of the invariant charaeter of the 
coefficients of an arbitrary R. Consider first the identities of the 


form 


(n-1) 


t 
Dio ..edel iti... ng (H) = Du) + RB. 


The expression D, (u) is of order n-l and, in fact, 


| 
Sow. ** e a 
ane eve oo , ¢ oa 
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Di ya 7 > “ y{n-2) Wee ve 
u) ' + wae te see. hee _ = 
i OR OK; Oy Ak, =? Jr jel DS oy, Dey hy a 
j#l 
Th 
+i (n- - 
1 Ke] ewe ee ‘K-) Ke] «> x ee + 
dei kel mia RPK eh 


j*k 


tere + AU, 


We may write this more canpactly as 


: = (n-1) = {n«2 ), = whe B+ coo t &2Ul, 
8) DC = Hay 2S 5) (43) F 2 iy 2, 23) (45k 1 
j#i J#1 K#i 
J#K 
To illuminate the notation, the superscript on u indicates the order of 
differentiation, the subscript on u, which is in parentheses, indicates 
the variables with respect to which u is not differentiated and the 
subscripts of the a coefficients, which are in parentheses, indicate 
the integer or integers auitted from the subscripts of the a coefficients 


in D. For example, 


ae) “2 66% Juk jtl..28 


(5x) “10 o°00 jel Jtlecc kel k4l..-n. 


Observe that a jk) =a hence the factor 2 for similar reasons, the 


(kj)? 
Moe term in this expression contains the factor 1 


q-1)§ 
If we scrutinize equation (8) we find that there is just one term 


of order n-l, and its coefficient is 1; there are n-l terms of order 
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n-2, and their coefficients are the a, 3) there are Cs) terms of 
order (n-3), with coefficients @(3,)3 and in general there are (San) 


terms of order meq. Consider now the effect of the operator 


(n=1) = 9 | (i) 
Pieesi-l it]... 0 Sz F O(g eee ete opps te By U).. 
4 | 


(m1) my = (n-1) ~2) 
wMasedel me. es 5 * 2 00" (3) ? (° oil * ea)" “w) "a" 
j#a OF 

m2 a 
(9) Gt 

e e ®(5)"(ge) * = 

th iced . ; 

jk 


! 


Equation (9) tells us immediately that there are no invariant coefficients 
for the terms of order n-l. This statement is readily seen to be true 
no matter which identity we consider. Thus a discussion of the invariants 
mist begin with the terms of order n-,. 

Fran the fact that the indices i, j, and k are distinct we con- 
clude that (43) %(4,) for any Jj. We further note that (9) contains 
(n=1) + (al) in-2) =(3 different mixed derivatives of order n-2, and 
these are all “the possible mixed derivatives of this order. Thus in order 


to obtain the identity 


“1 
ae 441 vee 2 Pg(t) = D(u) + R, 


i.e. to make the expression D(u) appear om the right hand side of (9), we 


mist (at the very least) add and substract to the right hand side of (9) 


the terms a 4 3)" tray This results in the invariant coefficient of ardsr 


Nad s 





BF 


‘ : ~q ve, & i,j =1,2 eeell ij. 
(10) Pty) oh 2) MET S(aH) ? aa 
x 
1 
Since there are n choices for i, and n-l choices for j, we have 
a total of n(n-1) invariant coefficients for the terms of order n-2. 
We will now prove the following: 
Theorem V: The invariant coefficients for the terms Ue 43) of order 


of any identity are of the form 





n-2 which appear in the remainier_ 


3265) +9(5)9(4) 7 254) ° 


Ox 


a 


These "invariants are ali true invariants, i.e. invariant under the 


ange of variables u= Na, no invariant coefficient for terms of 








Consider any identity given by a sequence of operators, 


2) sp? pO! pls) . . 
tase te “Iti e ee alas 5 aos” D(u) + R 


where (n-a)+(n-b)+... + (n-s) 
If we write this out in more detail, we find that such an identity is 





actually 
ee -a-2 
dg a, ed 122 ogo 9 
OX, OK IKQ jet xy: - -IKy---OXy, 2 fe x, OR OX, 
(i# lplz,---y4g) (i,jeipt,--, 1a; 1+ J) (4, jpK# 114° -yLas 14h.k) 
(n-b) n (n-b-1) 
ae ie s-fo 
| xy Oxy okn fi Md Je dey oe, |” 
(i# jys--+ode) uf # it : dee i+}) 
hai S) y je .. »| =eDOneeae 
E oy, 0%, + > ai) 5 ‘ os DKye OR, s 


(14 1,,--- 51s) oe 2 la 3 ied) 
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By an inspection of the left hand side of (11) we can determine how each 


term in D(u) is made to appear on the right hand side of (11). To 


obtain the ni order term On » We mist canbine the highest order 
AX, °° 0x, 
differentiation in each operator 38) pi?) pia) » and apply 
ints es tk Pie. J Ky oeoky 


) (n -S) u 


IX} 000 AX] coe OXy 
(i% L,-0051,) 


the result to ‘ (he determination of this n@4 order 


term is, in fact, the criteria used to determine which combination of 
operators yields an identity.) The terms of order n-1 which arise are 
obtained in two ways. They are obtained by combining the highest order 
(a) (b) vp @) 
differentiation in each operator Di eoedgs Djyses dy SP, Icy 280g? and apply- 


ing the result to each term of the sum 


22. 
in "(3) a Se, 
(22) $i semaEree 


(i,j =1,...2,3 1# 3) 


Using Leibnitz's rule we see that, among others, the terms 
nol 
(3) 805) J = 1,2, 000)m 
J . 1y..04, 


mist appear. Such terms of ordsr (m1) also arise when we combine, one 


by one, a differentiation of second highest order in any one operator 


pe) ft ; pia) y. With the highest order differentiation in every other 
000 cee oo8 q 


1 1 
operator pi?) 4 a 8aay Ke DE sor and apply this result to 


a 
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2d (n-5) 
(33) KI 
i#1,...1, 


The coefficients of the terms of order n-2, which are the terms 
we are concerned with in this theorem, arise in four ways. We may ob- 
tain such terms by combining the highest order differentiation in each 


operator Se Aly a eN and applying this result to each term 


of the sum (12). Again using Leibnitz's rule, we see that, among others, 


the terms 


O &( 45) ou (n-2) i,J = Lye, 00 oS i ¥ a 


Sx (43) 
= L,J#ljeoe], 


mist appear. If we combine, one by one a differentiation of second 


highest order in any ome operator p\2) p{a) » with the highest 
poredgscees ky ooek, 


f 
order differentiation in every other operator D{®) oh ie | and 
las a? ! Raye 3ck, 


apply this result to each term of the sum (12), we see that, among others, 


the terms 
(n-2) 7 
&( 4) ®(4) 444) i, J = 1,2, ooeyD} i ¥ j 


1,3 #101, 


will arise. Ina similar manner, if we eanbine a differentiation of third 


highest order in any one operator oS + <agee PE ce with the highest 


(2). (a) « and apply 
wey 1 = q 
this result to (13), we obtain the terms of the form 
: (n-2) 
“(j) “(ge) ? 


while if we combine the highest order differentiation of each operator 


order differentiation in every other operator D 
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p*) * pv, $ and apply this result to each term of the sum 
le ae o@e@°0 Seem 
s a 
> (5) 5s 5 — 5, 
e086 eo@ ) 
j,k=1 gs *n 
j #xk (1,5, #1) +0013 1#),k) 


we must again obtain terms of the form 
(n-2) 

ee GS). 
Comparing the coefficients of these terms of order (n-<) and recalling 
the procedure necessary to make D(u) appear on the right hand side of 
(9), we see immediately that the invariant coefficients for the terns 
_’, which appear in the remainder R must be of the prescribed form 
(10). 

Knowing that these invariants mst be of this form, we proceed to 
show that they are true invariants. We consider the equation D(u)=0 
and the change of variables wu =Aq » where \= \ (X) is not zero. Using 
the chain rule for differentiation of a product, and employing the notation 


of (8) we have 


n, 
AS & ia! . $ )» ae), 1 a 


Ni 
> — (n-2) 
iz, paz 4 


ax. ix *(43) Ty xs ' 
(14) od 
a). > ¥ 2) ga) 


br —(n=1) ) 
“22 W453) th 2 ¢ 
k 


a Gs) 


OOes 


U 
el Ox,o0%, (ijk 
j#4 ji kvl yo 
dé 1" 1,2, e00)0 
afP2)_ gt, - (2-3) 
+ s609259 = 1,2 coe th i# j 
(45) =i "(43) +2. Sh "(45x ic) I~? I-39 av 3 
k#J 


k#} 
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and similar other expressions involving derivatives of order less than 
(n-2). Substituting the relations (14) into the equation D(u) = 0, and 
dividing through by d » we obtain 


D(a) = + > (ar, y* son 


ann” + 


+e7,, dier\yar dod ) eH ..20, 
ta & (ey (4) ee fA) W043) 
¥j 


Let us indicate these new coefficients with the obvious notation 
a = + di r 
Pata) =) 7 
8 (43) > 243) T (4) se + geak. 1 oA : 
dh 0X4 0X5 


Then the new invariant coefficients (associated with D), for the terms 
u are 
(a5)’ a 
a - 2 
ae i)* 85) 84) 7 2¢a3) 9 
= 
and in terms of the original coefficients, these becane 


yi oe +12) Sageh)( 0) + ee) ; 
a4 “3 
fant Bebo ashy a TA) 


= (2%(4) +°(4) (3) 7 “(sa)h a iihalen “ial 
(9 
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To show that no invariant coefficient for terms of lesser order is a true 
invariant, we must first investigate further the character of these other 


invariants, 


E. We will now give a general method for the determination of the ine 
variant coefficient of any order term in an arbitrary remainder R, 


Consider any identity given by a ssaquence of operators. 


(a) _(b) (a) p(s) = D(u) + 8, 
bait Site j,°etp a oie 1, eee aa (v) 


where (n-a) + (n-b)+...+¢n=s) =n. 


We will denote for brevity 


on) pla =Ho)e 


Tyecod, Kj cook, 
Suppose we wish to find the invariant coefficient for the derivative 
a!) » When ¢,U yi?) is a term in the expression of) sgt 
ao, By 32% 


a each differentiation : » L1=1,.0m, never appears in more than 
one operator on the left i “diae of (15), we see that the left hand 


side of (15) contains the term 


(p) 
£3 (e,) Vex : 


and no other term on the left hand side of (15) contains al?) n 2 
ae? Se? 


a factor. Thus, recalling the procedure necessary to make D(u) appear 


on the right hand side of (15), we see that the invariant coefficient 
(p) 


for wv must be 


My « o alt, 


(16) AN 0g) ~ Any 5 ethy 
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Now suppose we wish to find the invariant coefficient far the 
(p) 


derivative “in, o« om,’ when there is no term in the expression (3) l (u) 
3 
which contains uP) nm 28 @ factor, We may assume that Dt v] (u). 
are p 8 
mb ee (p=) 


contains a term Co Uy, oe oly 4 for seme j=1,2,...p, where we set 


te 


(eo) 
fees 


u. This statement is chviously true for, at worst, the expression 
1 (u) must contain a term ¢,u, (although e, my vanish identically). 


We may then apply LeRoux's rule, extended to an operator with mn inde- 


pendent variables, to the product _¢ au? Pan e. Since the original 


“J 
operator D never contains any differentiation Q multiply, for any 
x 


HD eras » RB, the same is true of the operator {J e Thus our extension of 


LeRoux's rule for this type of operator is 


n 
(17) Diuv)=uDdl(v)4 > 9B D' (w)H1 a as ee 
a ay 2 BL Pn Os 49% “ak 


J#k 
and hence 
086 ian, ¥" “Alte Mealy $30" te) ee Myajee 
p-j+k 
(3) (p) 
fet + ot Mp j41 222%) Oty Ht omy Pot bamela 
(p) 


From this we easily conclude that the invariant coefficient for 


oe oll, 


must be 


(19) a 7 (e,) -a ° 
aa : 
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We must point out that formula (19) does not give the complete 


picture. For if fs) (3) is the identity operator, then it must be 


that Pues. (p) t hand si 
a6 Co=a » an at a appears on the lef side 
‘ M0 0 oll My oo oll, = 


My «ell, 
of (15). That is to say, there is no need to add and substract 
(p) 
a u 
TM} oo eM, me eek, 
containing a at in R. In other words, there will be no invariant 


to the right hand side of (15), henee there is no term 


ecefficient for iP ) « The same conlusion holds whensver 
<2, 
a 
a appears ag a coefficient in any ome of the operators D 3 
My «e oll, i,eoed, 


4 


pr) eee, pia) « This is apparent from the left hand side equation 
J coodh> ky oo eK, 


(4), amd the discussion which follows it regarding the determination of each 
term in D(u) which appears on the right hand side of (11). To summarize, 


whenever a appears as a coefficient in any operator pl) ’ 
My eoeM, Lj oo. a 


p\?) 9 e086 pis) » which is on the left hand side of (15), there 
Jzovedy Lier]. , 
will be mo invariant coefficient for ylP appearing in R ao the 


D 
right hand side of (15). Otherwise, there will be an invariant coefficient 


(which may be zero) appearing in R om the right hand side of (15), which 


will be given by formula (16) or forma (19) as appropriate. 


F. As yet we have said nothing about the character or number of quasi- 
invariants, i.e. the invariant coefficients for derivatives of u of order 
leas than n-2. As we-have pointed out, however, there will be no invariant 
coefficient for a derivative i? when the.coefficient appears 
| gee 9 ers 

on the left hand side of any such identity (15). From this fact we may 


easily compute the maximum number of quasi-invariamtsa associated with any 





9) 


identity (15). 

The total number of terms in D(u) which are of order less than 
nee, provided no coefficient is zero, is simply = (7) > if we sub- 
tract from this the total mumber of coefficients of derivatives of order 
less than m-2 in D(u), which appear om the left hand side of (15), 
assuming all are not zero, we will-have the mximm number of quasi- 
invariants,which can appear in R, on the right hand side of (15). For 
if such a coefficient would appear om the left hand sids of (15) » were 
it not zero, then, regardless of the fact that it is zero, there will be 


no quasi~inveriant corresponding to this term. This mumber, then, is 


(20) s (3\ - a (s+ CP) ret Pal 


1*3 


(n-a)} (n=b)+ ..0+(m-s) = mn, where if c¢<3, we set ¢ A¢ KO. 
£3 
For n 2, the number of quasi-invariants is obviotdsly equal to 


zero. For n=e3, we observe that an unique quasi~invariant appears as the 
coefficient of wu in each identity, hence there are a total of 12 quasi- 
inyariants for the entire third order case. 

For n>3, however, the picture becomes somewhat clouded. This 
occurs because while no quasi-invariant can appear twice in the same 
identity, the same quasi-invariant can and does appear in more than one 
identity. Thus the determination of the total mumber of quesi-invariants 
for a given order n seems to be a rather tedious, lengthy, and unre- 
warding problem. To obtain an indication of the enormity of the task in- 
volved, however, it is quite simple to obtain upper and lower bounds on 


this number, which we shall designate M(n). 
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An inmediate lower bound is a consequemece of the fact that an 
unique quasi-invariant results in each identity as the cosfficient of u, 
Hence N(m) is a lower boumd for M(n). Yo obtain am upper bound, we 
merely observe that the total mumber of quasi-invariemts associated with 
D must be less than the mumber which would exist if every quasi-invariant 
appeared in omly ome idemtity. Now with every ordered integral partition 


of Js there are associated 





m™ {%ap% ap ++ % 


B\f J \f J-XK, \[d(%, +a, ay 
ING) x 2 XK 3 i x. identities of the form’ 


(21) vd) pm Xa), (m-06 2) | pl™®= Kudu) = d(u) + R, (see 


equations(6) ana (7 )o 


With each of these idemtities are associated a maximm number of quasi- 


invariants, given by (20), namely 


n nej k & 
e@ £()-(5 ()+2 2 ay 
Thus if we sum over A, the set of all ordered integral partitions 
of Jj, amd them sum over j, we will obtain the number of quasi-invariants 
which would exist if every quasi-invariant appeared in only ome idembity. 
This number is 
£ n=l 2 
- Woy =D vy () 
Jal 


where 





(24) 0° (a)(ag') of 


Then for n 2 3, we have 





N(n) & Mn) & M*(n). 


Computing a few values for M(n) we find the following: 


n N(n)  M (n) 
3 12 12 

y 7h 362 

a 540 8510 , 


G, Now we are in a position to complete the proof of Theorem V. We must 
show that no invariant coefficient for terms of order less than n-e can 
in general be a true invariant. We shall prove this rigorously oly for 
identities in woth = 2+ @(,.)7 and shall only indicate the proof 
for other identities. 

Consider the invariant coefficient for u » Where p<a-e, 


p My © ell, 
is contained 


and m r for any 1=—1,2,..06¢6,p. Then c,u 
1 ca <2 see fzo My oe oly ay 


in D (u) and in fact 


Q 
il 


° “mee ,, x ; 
(25) en mT, 
ooo My 


“p = By 


where of course r must be properly ordered with Mi, g008 rly 
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Now comsider the change of variables u = d a, where A= d(%) 2 Oo 


We extend the equations (14) to cbtain 


- 


ee me = * 
OX +00 9%, Oy tte, m, °¢ el, sce 
> 7 e: a _ 
tes, am 





ptl “paiel pHi4l 
OP y ¥ 
a = + ceo 
O7n, 20. 27m, ASy ay . 


If we subatitute (26) imto D(u) = 0, and divide through by » 7 


wo obtain D(H) o The coefficients in which we are interested are then 


mn 
er Pn t>_ “nh om 1 dA 
ra ™P =a  D 4sptl 1 Pd a any 


A 
(27) we ' 
+ % a aC ) + eae vat aoo 
jy kapel vd “ped PN “a, ae a 


“~~ 
jt 
we * S ate Roce "A i 
j,k, | = pal oa’ Segoe d OX, DX, Ox, 
‘jacal 43 ae 
n=p 
tooog 1 3 IN ° 


3 


Aon a + OF, 





ae 


Reytomeny +B omar GoM 


jg=Pptl 
J#O 
ri 
+, . dd 
—. Ty «ol M4? i tee et 
ieee NO 
jek 


) N-(p+1) 


+ = where x =x ; 
0 S —- 2 d=. on ; 2 3 
d “mit OF 19 “ity ee “mt, ne 7 


and a = a 1 od. 
r fee 
Bra +S fe 
In all the above we must comment thet the subscripts involved must all 
be properly ordered in accordance with the conditions of (5). 


The quasi-invariant for Vans » «ot in this particular identity is 
D 


simply a= c,) < Omy + « oI, », oF 


Q= 09 (a ) + a y(a ) - 
— = 5 x (r) ‘“*, ~ eee 
ox, ary my » my me 
fhen the transformed quasi-invariant, for ca m? will de 
aa °" 
(28) Q= 9 (8 ts (a peak 
x, mt eer (x) 2 71°" Z 


Substituting relations (27) , we find, after a sanewhat tedious camputation 


that 





LOO 





Q=Qy+ 

Pa (2 Ty oe mmr) + 2(x)° I oe +m me “Sm, oo oll a) on + 

Jepte 3 4 

2 
SS (Ses mmr) + (,)%m oo oll, = mmm ) a 0 Xr + 
=e " vm y 2,0 Ta 
jek 
J»keP te 
(29) 1 ‘ 
n- (p+2) 

- (2 Gy + “(r)*(35) = *(ar)? oo 3 d , ; = 
7040 "S41. jt May Dyed Met 


We observe that the coefficient of each derivative of r in this 
expression is itself a quasi-invariant which arises fram the particular 


identity we have chosen. 


Then using an obvious notation we write (29) as 


Pe ee ee eee 
d Mn, NO%ny 9 ny 
(30) 
. (n-(p42) 
+2 8p $ 2 ore eee +2 32) a a 


\ 


Now if p=n-2, then Q = Q + Qn 


Yih 
FS 
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But QO, ( 
0H, 


on! t Se)" 1. A Geir) cate, 


=. (1) + 20x) ol - 479) O, hence 


r 


Q=Q as was proved previously in the more general case. 
If p>n-2, however, then Q=Q +a linear combination of other invar- 
iants arising from the same identity. Thus QQ in general, for the 


particuler identities fon ae D. (u)@D(u) + R. 
me Y-] me) 2 


Although we shall not prove it rigorously here, it easy to see that 


for p>n-e, no matter which identity we choose, terms involving a) ; 


“my 
Q* 
ones y cooy Will arises in the expression for any transformed quasi- 


3%, 0%, 

invariant, since these factors appear in the transformed equation D(i) =O, 
Thus, we assert that for any quasi-invariant Q, Q %Q, and in fact we may 
expect that Q=Q + a linear cambination of other invariants arising from 
the same identity. This is the justification for the name “quasi-invariant, 


and this canpletes the proof of Theorem V. 


H. Im order to see the problems involved in applying the cascade method 
for the n®2 order operatars we have been considering, let us consider a 
particular kind of identity once again. We will study identities of the 


form 


t _(m-1) oe ea 
(31) 2. Di eopel Penal) S (2 +9(,)(0)) = 
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2 (3 eye)" (3)"*¢ae) ) "Gey +22 a REY te ees 
(321) im “’ 


where Q is the quasi-invariant associated with wu 
1, ° ee i 1, e¢ i a3 


int Ds Pp =1l, eee 30 


(m=-1) 
Lecoetlr41...R 


Denote D (u) = O uga vu=ui ot 
= 


J 
= O%(r) +2/,.) 2(4)°" 2( yy) 2, and observing that D(u)=0, we write 
= ue (x) “(J) jn), ; 


s 


‘= 2 ) 
(32) ee my 


‘ua > + Fa 
D(a ) = a T3*(5r) Ay seed gg dyeeet yg : 
jar 

It would be most desirable if all the T's amd all the Q's in (32) were 
identically zero, for if this were true we would have succeeded in reducing 
the order of the original equation by one, resulting in a single first 
order equation, coupled with an equation of order (mel). Ideally by this 
method we would like to reduce the n’4 order equation to n first order 
equations, but this happy result will, umfortumately, occur very seldan. 
In the more probable event that mot all the T's amd Q's vanish 
identically, we would then hope at least to be able to cascade the equa- 
tions in the manner of the preceding chapters, until such time as the T's 
and the Q's may all be identically zero. To this end it will obviously 
be necessary to impose certain rather rigid conditions upon the coefficients 


of D(u). Let the hypotheses be 
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(a) T= . ool), AT = ine ne hae * cae 


z ra “rel 


=rpe#o, 
rie 3 
(b) all coefficients are functions of x, only, 


(c) all coefficients with an "r" subscript are identical. 


Then 
D(u')=2 4 : 
(33) ia ) PS (iit » ifr 


Solving u = Ou ir » Tor u, we find 


x 
xr 


Uu= oe Je (60°) rs Ox, t F(x) ++ -%_7, may eom)f 
d 
hence - (2, dx a .\ Ax 
“slat ) §®(r) = 
m s o% . . “T fo , w ax, + a 


i ee oF oe cen deien 
0 14° 3 142 


Seer, q a, ,oX, 
a ae rf (2 7 ves, + : 


where q is the number of non-zero i 


Substituting these expressions into (33) we obtain 


A “a D (at) aa) = at fo deen vas, 47 : 
Meo. OAL, l 
oe 


We then take the partial derivative of both sides of this expression with 
respect to Xn» to obtain 
B,_\dx 
§ (xr) P et i Q 
=. a(y)Pr(u ) + d (eG ye. One D(u')|= 
T } <r x 


a Je) ae x, oe ’ ipfr ° 
ian 2 





10h, 


Cancelling out the exponentials and collecting terms we find that we have 
an expression 
» (ane) zo Ovag 


th 


where ,D is ann” order linear operator of the form (S), but with 


coefficients differing from D. These coefficients we may compute to be 


a!) = a(r)- O log? 





° 
oe ’ 

a! 2 &0o° — as oe 3 whnere 1 =r for some p =; , eels 

12 Ly as n Pp 
a’ a, ; r ’ , a 

A, ecel 580 Ly eoel as 1 Pr. ai T 3 where 

=2 — n=2 + Le ee Ne? (r) 
2 


ifr for ATy p=L, e009 Tim2 6 


If now we compute the new invariants we find 


/ : 
x ea - = ‘¢ *) j a § 2 sas Poke rtly seen 
J ‘are c 
and Q!, . eT ga. nese 
t7 ee ely W3 ee ee o 
ES 


But by hypothesis (c), a5) r for all of r and all permita- 


| 1] oeeinn3 
tions 17 eed _3 ¢ web it yr for any p=ls eocy Nede 

Thus if the new invariants are all zero, we may reduce the order of the 
equation D(a") = 0 by one. If the new invariants are not all zero, the 
hypotheses of this method are again satisfied, and we may cascade the 


equations in the hope that the invariants will eventually become all zero. 


This can only happen if T=n o a( 3) » where on is some positive integer. 
Ox 
r 
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Other methods can be developed for identities of this form and for 
identities of other forms, but in every case, severe restrictions must be 
placed on the coefficients before cascading may begin. We see then that 
the problem of attaeking these equations in this manner becomes extremely 
difficult and restrictive as n increases. For n>h the problem becomes 
enormous, with staggering — of invariants and identities to be manipul- 
ated, and a large number of restrictions on the coefficients necessary to 
cascade the equations. If it is desired to cascade the equations it is 


recommended that only first order, substitutions be employed. 





SECTION VI 


i 
X 


‘ 


A BRIEF SUMMARY OF EARLIER EXTENSIONS 


A. Extension of the Laplace cascade method have been made by many note 
able mathematiciams. One of the earliest of such extensions was made 
by Darboux himself > Darboux considers the following system of equa- 
tions of second order, with m independent variables. 

Let Po J Q 1 O. go’ Sin be a system of n independent variables. 


and consider the system of equations ; 
(1) 20 x, = 4, ou +g s 15k=0,1,2,..=-1; ifk, 

Px 
This system contains n(m-1) limear hyperbolic equations of second order, 
and we observe that this Z an overdetermined system, whereas the system 
Cedsiasrea by us in Section III is exactly determined. Darboux seeks 
first a mecessary condition to give n linearly independent solutions 
in addition to the trivial solution u= constant. To accanplish this, 


Darboux forms the third derivative by taking & of both sides of (1). 


Interchanging the indices | and k, and equating coefficients of like 


derivatives of u, we obtain the relations 


(2) = =» Se ¢ tale gtk (l=k-}). 
Interchanging tnd indices i and | does not change the right hand side 
of (2) and hence 

da4~ - O%K 

20, ay 
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Holding k constant, we see from integrability conditions that there 


must exist a function, call it logH» such that 


(3) ai, 1 OF (iv ic) | 


H,. oPi 7 


Then condition (2) becomes 


(4) om. 1 OR OR ae : 
300, aq, dei #8 30, 


and the system (1) bas the form 





(5) ou 2 MQ. OH ibn oa 
2028, Mc OG; OQ, Fi dex Ber : 


1,k=0,1,2, eo oy Rl=L3 i¥k, 


Let fy (Wed u ~ 1 om du - 1 O81 du 


age, Bk oes 1 Ok Hy den 304 


Using condition (4) it is easy to verify that 


(6) i (H,)=0 > iE L, 


To start the Laplace method we assume that we have an equation 
of the form (5), satisfying the conditions of integrability (4). Thus 


we consider the particular equation (for fixed i, k), 


(7) £4y(4) mee 
Derboux then considers the substitution function V defined by the 


relation 
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OH k(u +V). 


(8) du = 1 O08 
By 


i 


If we substitute (8) into (7), the equation then becanes 


10%, gv 4[d Loam. 1 JE DH | ¢ _ 
(9) Re Pr Oy * PPey, i See 3p. (ut¥ )= 0, 


We see that cur invariant hy is now 


“nr ( ~ 1 08, a a, 
wae oe 3 


and if h,,=0, the equation (9) has the simple solution 


V = constant , 
We may gain a little insight into these proceedings if we adopt 


the notation of Section V. Thus, for equation (7), 


Da 3° 1 A o - i De 
Ce i 1 Px Br Mx Ay 


d tm 2 
ee Oi ik ¥ e's 3B; 
e, 


Then the substitution function V is defined by 


(10) ae (u) - ,°V.> © , 


while the invariant h,, may be written 
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__ y 
(11) +f. Fen. 0, 30? ; 


amd hence (9) becomes 


/ 
D 


(12) tk ( *1)] u é 


/ 
D (a) = D,,(u)+(D 
bi in 1k ( ato, 


If -a 7 =0, we see that (10) reduces the order of 


: ( 
D 
ik 
G 
(12) by ome, and the equation can be solved by quadratures. 


In the more likely event that ig Fs) # O, we may cascade the 
‘ i 


% 


equations in the following manner. Let 


b= = Sk  Dy(By) 
a 


Oy, 
oH, °* 


eT 


- 


Then 


k H 


ov. & (us) | 


k 
and V satisfies the equation 


(13) ov - 1 We WV - 1 Wi w= 
80:90 he 00; A, tO, 9D} 


| | = 
In like manner, if we introduce the quantities L' defined by the 


relations 
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bre B, Ox! - a! OF ; (x #4 #k) , 





then the function | , for distinct values of i' and k', satisfies 


the system of equations 


(14) + a 20% ww yi Okt DV 

Ox'9Pxt Ty" D@y' Oa" Txt O04" Py? 
which is of the same formas (5). -Hence we may iterate the substituion, 
and cascade the squations until such time as the corresponding invariants 


may vanish. 


22 
B. In 1899, J. Le Roux! extended the method of Laplace to linear 


partial differential equations of greater order than second. Le Roux, 
however considered the equation with only two independent variables, 
while we, in Section VY, have considered the equation in n independent 
variables. Im the notation of Le Roux, the equation to be considered is 


the nrh order equation 


X+B 
(aid) D(z) = 2. n | hn a2 z=, 
«1 @l(n-a-e)! ** Ax%XDy® 
where 04 6 En, XF n, 6# nu. 
If we supposs that the highest order of differentiation with respect to 
x im equation (15) is squal to m-p, then the collection of terms 
which contain such a differentiation are 


>) nepti ) n-p+i-1 


m=-P 
= + b ed a Oe + oo + & d a Z 
- ee +! ; jx*®) z= =~ J. ) n-p & 


x 





We denote the expression 


oO tro te 


Oy 1 yi a1 
g*# 
as the differential multiplier of the term zZ , and we further 


Ox? 


assume that the coefficient a is equal to ome. Le Roux cossiders the 
equation (15) from the point of view that there exists a particular 


integral of the form of Euler 


coe tude 


where X is an “arbitrary"function of x, and the coefficients u 


are functions of x andy. 


If we regard D(z) as a polymanial in O and 2 we may introduce 
the notatiom of Section V. We designete oa = aes expression ob- 
xy 
tained whem we differentiate D with respect to 0 p times, with 
respect to q times, and apply this operator ae variable Z. For 
y 
iastance, the (peg) ©” derivative of the term 
tea 
“ag OxOdy* 
is K-Pp+B-q 


x(«-1) ---(«- pti) 6(@-1)---(@-q+l) A “65 RP) yea 


and hence 


(p+q) -p+@-4 
(17) ®D (z)=n(n-1)...(m-p-q+1) > (n=p-4s A ar Zz 
"y (=p)! (B-g)n-w-e)! “QR EF, HA 


ay 


a 





ee 


Since (15) is asswusd. to have order (n-p) with respect to x, we see 


that the cosfficients of (16) must satisfy the equations 


(np) 
apy ae (to) = ° 
(18) (n-p) (n=p-1) 
yr = YY eau wa els 
(n=p) — 
1 9 n-p-1) me 
ene =” ik + apa mee (ig) AP <a 


@ 2 @ 
a 


In the case when the variable x is a simple characteristic variable, Le 


Roux considers the transformation — 


gaa 1 Dna (z) = O22 +nA Z9 
: nel)? x" Oy =o 
and defines the functions V and v, 
nA dy 
- 2 =29 n=l 
(19) 2 
19 zZ nA ay « 
= al i 4 nel - ov 
2, 1 7° 


We then consider the equation resulting fram (15) by means of the trans- 


formation (19), and collecting on the right all terms which do not contain 
any differentiation with respect to y, we obtain 


(20) AW,)=), ots i, nt + d a Sipe tothe 


where /\ (V, ) designates a Ce aaa expression of order (n-1), 
which the coefficient of a Was is equal to unity, The coefficients 


are given by the formula. 


mA dy (p=1) (p~i) 
(21) d= 24 en Ayr D pel (tp) = 3 ety DaPri Uo Cae), 
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The order p of the right hand side of (20) is, in gomeral, equal to 
(n-2). It will be less oly if 


Dinca (ug) = 0, 


Le Roux then proceeds to define the imvariemts h, by the follow- 


‘ing rule. If the A's regarded as fuactions of y are all linearly 


independent, then 


Xo Ai woe OH 

ddo Do. we dd 

0y ov Oy 
(22) n= opis anf (uh 


> 0.AW se 


which are (pl) in mumber, sotting h=),. If the d's, regarded as 
functions of y, are mot linearly independant, we define de, as the 
first K which is linearly independent of ro A, as the ‘faret 'Aahich 


‘ds linearly independent. of A,cana. Ay, etc. Then the determinant (22) 


is replaced by 
» ec, ed ro, 


aye dor ro; 


0 


(221) h, = Oy a oy 
hs rec 3 Docs 





dy’ oy dy 


To show that these are generalizations of the Darboux invariants, we must 
consider the transformation 


222. f(x,y) , 


(23) \" =  (x') 
| y= iY (x',y') c 
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The first of these transformations is effected without changing the value 





u 
of V , sipce u' = ° » and hence 
2 2(x,y) 
2! fae @ 
ut ‘ 
° eS 


The transformation of the independent variables, however, causes equation 


(20) to became 


P 
(24) A Wo= Xo JV +d ons set ApY 
f. 


Ox ox 
where 4 S : 
oy v= YU, se 


a XK eed, @x) 2 2 
oO | © Z ; are 
t === 
(25) Ar > (Ba) "D+ P19 0) 


Ox! 


Nady /m: \Giet-pt2) 
; (3) ; (hs So AtE,0 d,) 


the coefficients © being dependent only on x. If we use the relations 


(25) to compute the new invariants hj » we will find that 


(441) (i+2) (m-2-p) (141) 
ni « Q esa RIE 
7 (ge, 2x) a oe 24 


Since in gemeral pen-c, the above relatim reduces to 


i*])( i+2 
Mt = (2 2a. 
ay’ Oo 


Thus the determinant h, is reproduced mitliplied by the (i+1) (i+2) 
2 


a: 
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Power of the Jacobian of the transformation. This justifies the name 
"invariant" for the determinants h» and LeRoux calls these the first 
generalization of the Darboux invariants. 

Now for the equation (15) to admit an integral of the form 
Z@u~, X being an arbitrary function of x only, it is necessary 
and sufficient that the coefficients A of (20) are all zero. This is 


evident if we but consider 


tt tr oat (n-1) =1 
(ux) X(D(4.) + E'DL(a) + KE" Diol ug) + ++. hoes Dine (Yo)™ 


s 


Since X is arbitrary, and this equation mst be identically zero, it 


follows that 
t ef nel =— 
D(u,) = D,(u,) = D_2(Uo) Miuee = na14o) = w; 


which proves the assertion. Thus we see that the vanishing of the invar- 
iants plays the same role for (15) as it does for the second order 
equation. 

LeRoux then considers the case when x is a multiple charactoer- 


istic variable, so that (15) has the form 


Nii Nem —| 
D = P 3 O | = 
ee > aa Va S rat ot Ge LaF O , 


where dD ' d 1)°°‘' designate the differential multipliers, o.g. 


@.= Say + ar 1. Soe * "2 a San t eco ' (xr mu, Le 


The first substitution we have already considered, namely 


1 = % S(E), 
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We attempt to decompose @o(z) into differential factors, so 


that D(x) FX Hy oe Xray Uy O a. 20 Le 4 1 Oe 22 


al. 
OF Xe OF Kee Oy MDF Ty 
Then we define 


= 


Uo 


ps 
on 


<4 | 
x) Oy “1098 zhe y 
2 


0 
UW i = & ew & U 0 20 
a "V2" Fe OF dy X, Y WU ° BY HK ot, ) 


2 
oy 

5 =m; 2 L 
Q 


aad using these we are able to reduce the equation. If @ (2) does not 
have such a deccmposition, we may still reduce the equations by a method © 
which will be described briefly. The proof of the following is rather 
lengthy and detailed amd does mot bear repeating here. It may be found, 


however im LeRouwr's work previously noted (footnote 22). 
Let uy Ys Uoyere demote a set of limearly independent integrals 


of the equations (18) which define the conditions om the coefficieuts for 
the existence of a solution of the form of (16). We then define the 


functions 


0%, ou a 
U : Uu © @ ©. 
Pome | ° oF y2 © | ow oe | 
aL. ee U 0 My 3 2 u 3 ) ay 
=o ik © MH = Che am 
| Oo” . “Oxe 
We = oF Ys 





nts f 


Set 
1, m Kk, 5 


a | 6 
Kyo1 
The equation 


) 8 
| os wh oF ld aa 


admits as solutions the functions old | cools » Then the set of treasfor- 


mations is easily defined by setting 


yew 2 (=) 
2 %>li $i 39 (2) hi Z(t) 


a A 3 ° 3 O 
a Bu I, ee +7, 1 yy i 2 (&). 


If there exists a particular integral of the form of Euler, the invar- 
jants ] will beeane zero after a certain number of iterations, and 


the chain of transformations will stop. 


| : 2 
C. Following the work of Darboux and LeRoux, Dini ° made an extension 


a 
of the Laplace method to the limear second order equation im an arbitrary 
number of independent variables. Dini considered the equation 
n a 
2 
a A os G, Iz 4Net HHO, A,,= 
} - 
i OX 0, j=l J Ox tJ Jt, 


and posed the question, can we transform (27), imto the form 


(27) 


n 
a9 7 
k +MOt = Oz = 
(28) 2 i Tx, eZ 7 Se Br, 


where 
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“3 


and the functions a1, 8, +os,8,) Dy Ky,e0ekg, M are (2n+2) functions 


Zl 
(29) OF 2. 4 eye™ 


to be chosen arbitrarily? If we assume that this has been accamnplished, 


substitute (29) into (28), and equate coefficients with the left side of 


(27), we optain the (mtl)(nt2) equations 
2 


(30) kyS5 + j= 2A,, 5 1,JF1,..0,R , 
n 
(31) ky + aM4+5” ky dag +X, = G , Lyssa , 
j=l Ox 


J 


A. 
(32) Mt > x, db y+LaN 
gear (5 


for the (2n+2) umknowms. We observe, however, that im (30), the k, 
and a j appear only as products kya 3? and not all of them are zero if 
(27) is to be effectively of second order. Thus we may select. arbitrar-~ 
ily ome k, say kw and set k=, without affecting these relations; 

and hence there are actually only 2n-1l unkaowns a, and kK, If we consider 
first the equations (30), we note that for n=2, these are 3 equations 
in 3 unlmowns, and this system is determined. For n& 3, however, the 


system is overdetermined, and one would expect to have (n-1) (n-2) rel= 
2 


ations among the coefficients A,, so that only 2n-1 of the equations (30) 
will be independent. 

We ses then, that before we may even begin to seek apply a 
"Laplace method” to equation (27), the coefficients As j mist satisfy 


(m-1)(n-2) conditions. Dini showed that in general these conditions may 
2 
be stated as 





at 


2 2 s 
ee ne ae? S, Ce ane - Apher VAng - Anndss , "28 = or3y 000275 


Ir#4s 
J 
where € = £1,€,=21, ani these signs are determined by the sign of 


the radicals of the roots of the quadratic equation 


Ady - 2A... Bt tye . 
Here Ann is considered to be the first non-vanishing coefficient 
Agy: 121, ooeyRe If all the A,, are zero, the equations (33) are 
satisfied identically for all Ags, 1=1,..0,h. 
The classification of the various allowable types of equation 

(27) becomes the next important consideration. We list here without 
proof same of the more impavtant properties of (27) derived by Dini. 
First, if n>2, and all coefficients of (27) are real valued, ther 
for any pair of variables (x,,x,) for which the partial differential 
—- is not parabolic, the equation mist always be of the same type. 
That is, ignoring all pairs of variables for which the equation is para- 
bolic, and considering all remsining pairs of variables, the equation 
mist always be elliptic, or always be hyperbolic type, with respect toveach 
pair of variables. Another property, true even for conplex-valued coef- 
ficients: if Ay), #0, and the partial differential equation (27) is 
parabolic with respect to the pairs (x),5%,.) and. (x,2%,) » then it will 
also be parabolic with respect to the pair (x,9%,,) © These properties 
follow from the conditions (33).- 

| Dini then defines the following terminology. Equation (27) is 
said to be of parabolic type provided that it is parbolie with respect to 


every pair of variables (xp2X,) 53 if there exists at least one pair of 
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variables (x,.9X,) such thet (27) is not parabolic with respect to 
(x,,X,) then (27) is said to be of elliptié (hyperbolic) type if it is 
elliptic (hyperbolic) with respect to every pair of variables (x»)%g) 


arabolie type: otherwise (27) is said to be of 





mixed type. Then Dini proves that if (27) is not of parabolic type, 
then either there exist exactly two systems of values of k and a which 
satisfy (30), iatere exists none, while if (27) is of parabolic. type, 
there exist at most one system of values of k and a satisfying (30). 

In the first case, one obtains one, system of values from the other merely 
by imterchenging the k's and the a's. Denoting these values by (k,a) 
and (a,k) respectively, we may call them conjugate seta. 

Now our first uae im attempting to reduce (27) to the inte- 
gration of two partial differential equations of first order, is that 
perkaps (28) will reduce to a first order equation in © alone. Im ad- 
dition to the (m-1)(n-2) cosditions previously noted as necessary to 
redues (27) to: (58) and (29), this requires n-l further conditions 
to insure thet «,= Xx 2 coe X= L = 0. These conditions exise fram 
the (n+1) equations (31) end (32) which determine only the two un- 
knowns b and M. Assuming that the equation is not of parabolie type 
and hence, without loss of generality, that it is not parabolie with res- 
peet to the pair (x,%,) » Dini shows that the nel plitasciens ean be formulated 


as 


(34) N- Mb = BO 





ae 


aaa aon lee soa 
(35) ve ao ' Fs A, = 0 3 5 = fel coeyh 
A A Ge A 
1s 2B gi 8 
ma ya A 
where A_= a ky. 28 , and Be y kK ob : 
rel 0X, rol OX, 


We may look upon L and the OC : im (28) as playing the role 
of the invariants. Thus to say that the conditions (34) and (35) 
are satisfied, is to say that the invariants of (27) are all zero. 
If we denote D*(u) = M, + ae is j ou , these inveriants may be written 


j=l dX, 
as 


= ii p'(b), 
(36) 


Y= 1) aa 
i G,- k,b D (a,) 


Im order to start the chain of oquations when these imvariants are not 
all zero, Dini imposes conditions similar to those imposed by us om the 
third order equations of Section IV. For equation (27) to be cascaded, 


Dini requires that conditions (35) be satisfied. That is to say, the 
invariants XK, are all zero, but the invariant L is not. Then 


equation (28) may be solved for «2, im terms of © and its first 
partial derivatives. This expressiom is substituted into (29) anda 
meow second order equation, for © , is obtaimed. Not only will this 
equation be of the same formas (27), but the second order coefficients 


A, j are reproduced. Hemee our new equation is 





Lee 


n 2 n 
(37) 2 Ay, a 4 > G 08 +r 4 g's 0 
i,j Fl XOX, j=l Ox, . 


The new coefficients are given by 


Q 
G = G@ *# K - AwEK L 
8 8 6 .& 


8 8 J 


(38) y= N +> a Om - m, -2 
0% ; 


H = H(b +H, - La) 


k 
where K = 2. a 9 8, and H.? L, are given by 


n 
= Ol 
. >_ Saeed 4 
vex] x 
Sines the relations (30) used to determine (k,a) for (27) involve 


only the coefficients A it immediately follows that if there exists 


ana 
a set (k,a) which splits (27) into (28) and (29), the same set 
(k,a) will suffice to split (37) analogously. It will be necessary 
g 
however to compute a new M and b , and a new set of invariants L , 


g g 
and eC; . The conditions under which the % , are zero becom 


ou aan @ - Ay 
g 
G. « 
(39) Alo ‘50 2-85 Eo , for pe3eh cso 
A A G . A 


13 2a 8 8 
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which are quite similar to the conditions (35). If these conditions 
are satisfied, and further Le N -Mb._ - Bee) then (37) reduces 
to two first order eee as desired. If L' is not sero, however 
we must again require that eqhiitt ides (39) be satisfied before we can 
iterate the process to obtain a third second order squation of the form 
eile 

Dini's very remarkable result, however, is the following. If the 
conditions (35) and the conditions (39) are all satisfied, then we 
need impose no further condtions to iterate the process indefinitely. 


That is to say, if the invariants %, and 4 ; are all zero, then all 


invariants oc, @) will be zero for any positive number m of iteratians. 
This result becomes evident if we compute, say, G. by the rule given 
Hl ( 8 9 4 tf A ft 
n (38) for Gy substitute G, -A4A,, Gp 


corresponding terms in (39), and employ the rule for evaluating a deter- 


v9 
A,, and G -A_ for the 
s g 


minant when one of the columms is a sum of two or more columns. Each of 
the resulting determinants in the sum must vanish as a direct consequence 
of conditions (35) and (39). 

In summary then, in order for the equation (27) to be cascaded 
as many times as necessary for the L invariant to vanish, there must be 


satisfied a total of 


(n-1)(n-2) + 2 (n-2) = (n+3)(n-2) 
2 2] 


conditions on the coefficients, The first (n-1)(n-2) conditions are 
2 
necessary to reduce (27) to two first order equations, while the remain- 


ing 2(n-2) are necessary to permit the cascading to continues. We observe 


that for n=2, no conditions are necessary, and this is the original case 
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considered by Laplace, Legendre and Darboux. 


(2h) 
D. Finally, mention should be mado of the work of Burgatti who 


considered the equation of elliptic type, 


Doe a caf 1h Okc Be O, 
(40) 52°93," sat —e 


Burgatti found that the two expressions 


by 
i 
Qy 
> 

5 
QJ 
ov’ 


Y 


(41) 


4. ar + be a 2e 
eo 


fe 
tt 
le 
is 
le 


are invariant relative to the transformation z =) z'. Further, if 
H and K are both zero, equation (40) may be reduced to 
dz : eo. 
+ GE a0, 
Pe Oy 


which is Laplace's equation; if H is zero but K is not zero, then 


t 





(40) takes the form 
QPa' 4 Dez te'z'=0 s 
axe s° Oy y" 


while if K is zero, but H is not zero, (40) becomes 
2 t 8 ') 
ox ay 


where O¢ and @ are functions of x and J. 


In the notation of Section V, let 


pz O° 49° 409 4nd +e, 


oe 


Ox dy? 9X OF 





12) 


then mh= 20 + 2 
Ox 

DD = 29 + hb 
oy 


and the invariants (41) then become 


| 
H = Dye) =D, (>) = (Dye) - ©) - (Dg(d) = 9) 


(42) a ]. 


jw 


/ 
kK = p/(a) dp D,(b) - C= (D(a) -¢c)¢ (D,() a) 
2 2 2 











Since the change of faye and the vanishing of invariants does 
not reduce the equation to two first order equations, no consideration 
can be given to cascading the equations of elliptic type utilizing the 
invariants H and K. The work of Dini shows us, however, that this 
equation can be cascaded if we admit the possibility of camplex coeffi 


cients. 
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